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16749 BLANEY-CRIDDLE COEFFICIENTS FOR TURF GRASSES 


KEY WORDS: Evapotranspiration; Irrigation; Lysimeters; Turf grasses; 
Water consumption; Water supply 


ABSTRACT: Evapotranspiration rates were measured for several turf grasses 
Arizona, California, Colorado, and Wyoming. These data were used calibrate the 
Blaney-Criddle method for estimating evapotranspiration rates. Warm-season grasses 
used less water than did cool-season grasses when grown warm climates. However, 
Kentucky bluegrass did exhibit heat stress that reduced its evapotranspiration rates 
below what would expected when soil temperatures exceed 25/DEG The data 
were collected urban areas and should representative city-wide water 
requirements for turfed areas. 


REFERENCE: Borrelli, John (Assoc. Prof. Agricultural Engrg., P.O. Box 3354 
Univ. Wyoming, Laramie, Wyo. 82071), Pochop, Larry O., Kneebone, William R., 
Pepper, L., Danielson, Robert E., Hart, William E., and Youngner, Victor B., 
“Blaney-Criddle Coefficients for Western Turf Grasses,” Journal the Irrigation and 
Drainage Division, ASCE, Vol. 107, No. IR4, Proc. Paper 16749, December, 1981, 
pp. 333-341 


16740 SOIL PERCOLATION RATES 


KEY WORDS: Field tests; Infiltration (water); Irrigation systems; Nuclear 
probes; Percolation tests; Soil moisture; Subsurface investigations 


ABSTRACT: study was undertaken investigate the feasibility using subsurface 
neutron moisture probes estimate soil percolation rates. The measurement procedure 
developed based upon the subsurface soil moisture content versus time relationship 
obtained with subsurface neutron probe soil column which one-dimensional 
percolation occurring. The results from ten laboratory percolation tests employing 


column well-graded sand were used validate the measurement procedure. The 
proposed procedure appears acceptable for homogeneous deposits well-graded 
sand. The validity the procedure for layered soil profiles shallow groundwater 
table conditions must investigated. The procedure proposed simpler than 
traditional methods (e.g., lysimetric) measure percolation rates. 


REFERENCE: Vinson, Ted (Assoc. Prof. Civ. Engrg., Oregon State Univ., 
Corvallis, Oreg. 97331), and Mahar, Larry, Percolation Rates From Neutron 
Probe Measurements,” the Irrigation and Drainage Division, ASCE, Vol. 
107, No. IR4, Proc. Paper 16740, December, 1981, pp. 343-359 


16747 DIMENSIONLESS STREAM ADVANCE SLOPING BORDERS 


KEY WORDS: Border irrigation; Graphs; Hydraulics; Mathematical models; 
Open channel flow; Surface irrigation; Water distribution 


ABSTRACT: The optimum choice characteristic reference variables used put the 
zero-inertia governing equations continuity and momentum with boundary 
condition, into dimensionless form not obvious. The effect different choices 
noted, are the effects choosing different formulas for field roughness and 
infiltration. The choice normal depth for characteristic dept, characteristic distance 
equal the quotient normal depth and bottom slope, and characteristic time equal 
the time travel the characteristic distance normal velocity leads useful 
two-parameter set dimensionless curves for advance prior cut off border 
indefinite length. These are presented for series Kostiakov-infiltration-formula 
dimensionless coefficients and exponents. proves possible present virtually all 
practical field and laboratory combinations input variables—inflow rate and border 
slope, Manning roughness, and infiltration—in ten graphs, each spanning log cycles. 


REFERENCE: Strelkoff, Theodor (Independent Consultant, Liberty St., San 
Francisco, Calif.), and Clemmens, Albert J., “Dimensionless Stream Advance 
Sloping Borders,” Journal the Irrigation and Drainage Division, ASCE, Vol. 107, 
No. IR4, Proc. Paper 16747, December, 1981, pp. 361-382 


U.S. Conversion Factors 


accordance with the October, 1970 action the ASCE Board Direction, which stated 
that all publications the Society should list all measurements both U.S. Customary and 
(International System) units, the following list contains conversion factors enable readers 
compute the unit values measurements. complete guide the system and its 
use has been published the American Society for Testing and Materials. Copies this 
publication (ASTM E-380) can purchased from ASCE price $3.00 each; orders must 


prepaid. 


All authors Journal papers are being asked prepare their papers this dual-unit format. 
provide preliminary assistance authors, the following list conversion factors and guides 
are recommended the ASCE Committee Metrication. 


convert 


inches (in.) 

feet (ft) 

yards (yd) 

miles (miles) 

square inches (sq in.) 
square feet (sq ft) 

square yards (sq yd) 

square miles (sq miles) 
acres (acre) 

cubic inches (cu in.) 

cubic feet (cu ft) 

cubic yards (cu yd) 

pounds mass 

tons (ton) mass 

pound force 

kilogram force (kgf) 

pounds per square foot (psf) 
pounds per square inch (psi) 
U.S. gallons (gal) 

acre-feet (acre-ft) 


Multiply 


millimeters (mm) 
meters (m) 

meters (m) 

kilometers (km) 

square millimeters 
square meters 
Square meters (m?) 
square kilometers 
hectares (ha) 

cubic millimeters 
cubic meters 

cubic meters 
kilograms (kg) 
kilograms (kg) 

newtons (N) 

newtons (N) 

pascals (Pa) 

kilopascals (kPa) 

liters (L) 

cubic meters 
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25.4 

0.305 

0.914 

1.61 

645 

0.093 

0.836 
2.59 

0.405 

16,400 
0.028 

0.765 

0.453 

907 

4.45 

9.81 

6.89 

3.79 

1,233 
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BLANEY-CRIDDLE COEFFICIENTS FOR WESTERN 
TurF GRASSES 


and Victor Youngner’ 


INTRODUCTION 


Studies the consumptive use water various crops are numerous. 
For example, Jensen (6) provides extensive review the status these 
studies. obvious from this review that nearly all consumptive use studies 
have dealt with agricultural crops. Very little information has been acquired 
concerning the water use turf grass, and information which does exist 
generally difficult locate (8,11,15). 

Turf areas are heavy consumers water and compete for water traditionally 
used for other purposes. Water applied turf obtained from municipal sources 
and thus given priority over other uses such agriculture and industry. 
Furthermore, most turf areas are thought being relatively small plots 
land, although there are over 8,100,000 turf the United States (2). 
these 8,100,000 ha, about 3,000,000 consist residential lawns (1,3,9). 

The consumptive use plants usually estimated through the use various 
empirical formulas. The formulas generally require some type local calibration 
order produce reliable results. Very few attempts calibrate these formulas 
for turf have been made (12,14). fact, Quackenbush and Phelan (12) state: 
research consumptive use lawn grasses urgently needed 
However, since their publication 1965, only the attempt Tovey, 
al. (14) calibrate consumptive use formulas has appeared national 

"Assoc. Prof. Agricultural Engrg., P.O. Box 3354, Univ. Wyoming, Laramie, 
Wyo. 82071. 

Agricultural Engrg., Univ. Wyoming, Laramie, Wyo. 

Plant Sciences, Univ. Arizona, Tucson, Ariz. 

Prof. Soils, Water, and Engrg., Univ. Arizona, Tucson, Ariz. 

Agronomy, Colorado State Univ., Ft. Collins, Colo. 

Assoc. Prof. Agricultural and Chemical Engrg., Colorado State Univ., 
Ft. Collins, Colo. 

Plant Sciences, Univ. California, Riverside, Calif. 

Note.—Discussion open until May 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication January 20, 
1981. This paper part the Journal the Irrigation and Drainage Division, Proceedings 
the American Society Civil Engineers, Vol. 107, No. IR4, December, 
1981. ISSN 0044-7978 /0004-0333 
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publications. Both these efforts considered grasses which were grown 
open agricultural fields rather than urban setting. 


Data Sources 


Measurements the consumptive use urban lawns throughout the western 
U.S. have been gathered and are used herein calibrate the SCS (13) version 


TABLE 1.—Site Descriptions 


Mean Mean 
annual annual 
tation, ure, frost 
degrees free 
Elevation, centi- centi- period, 
Location Latitude meters meters grade 
(1) (2) (3) (4) (5) (6) (7) 
Laramie, 
Wyo. 
Wheatland, 
Wyo. 
Fort Collins, 
Colo. 
Northglenn, 
Colo. 
Tucson, Ariz. 32°13’N 
Santa Ana, 
Calif. 


TABLE 2.—Summary Data Used 


Location Years data Grass type 
(1) (2) (3) 


Laramie, Wyo. 1976-1978 Kentucky bluegrass 

Wheatland, Wyo. 1976-1978 Kentucky bluegrass 

Fort Collins, Colo. 1977-1978 Kentucky bluegrass 

Northglenn, Colo. 1977-1978 Kentucky bluegrass 

Santa Ana, Calif. 1966-1969 Common bermuda, St. Augustine 

Santa Ana, Calif. 1970-1974 Kentucky bluegrass, tall fescue 

Tucson, Ariz. 1977-1979 Bermuda, Santa Ana and Tifgreen hybrid 
Bermudas, tall fescue, Zoysia, St. Augus- 
tine 


the Blaney-Criddle formula. Information describing the locations, years 
data, and type grasses given Tables and Consumptive use defined 
herein the evapotranspiration rate experienced when water availability 
never limiting factor. The grasses all cases were clipped turfs. 
Measurements Wyoming and Colorado were taken using small weighable 
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lysimeters located directly residential lawns (8). The lysimeters measured 
30.5 diam deep. They were made from polyvinyl chloride pipe 
and were surrounded retaining sleeves made from the same material. 
Consumptive use was measured weekly removing and weighing the lysimeters. 
Water was then added bring the soil near field capacity. Any excess 
water above field capacity was allowed drain storage reservoir gravel 


TABLE 3.—Mean Monthly Temperatures During Study Period, degrees centigrade 


Santa Ana 
1966- 
(4) (5) (6) (7) 


Month Wheatland Fort Collins Northglenn Santa Ana 
(2) (3) (4) (5) (6) (7) 


the bottom the lysimeters. assure that adequate moisture was always 
present, water was added during the week required that available moisture 
depletion was never allowed exceed 35%. The lysimeters were covered 
the homeowner during irrigation the surrounding lawn, thus preventing 
unwanted water from being applied. Detailed descriptions the construction 
and operation the lysimeters are given Pochop, al. (10). 


| 
Apr. 4.3 8.9 19.5 14.7 14.5 
May 8.4 13.3 22.5 17.5 16.6 
June 14.9 19.1 17.9 19.4 29.1 18.9 19.5 
July 18.6 22.8 22.6 23.7 30.8 21.9 
Aug. 19.9 20.1 21.1 29.6 23.6 22.4 
12.5 16.4 17.7 19.2 27.7 20.8 
Oct. 5.5 9.7 8.6 23.2 18.8 18.2 
TABLE 4.—Monthly Percentage Daytime Hours, Year 
Apr. 8.97 8.99 8.76 8.80 
May 10.09 10.13 9.62 9.71 
June 10.19 10.25 10.13 10.08 9.59 9.71 
July 10.31 10.36 10.26 10.21 9.77 9.88 
Aug. 9.60 9.62 9.57 9.55 9.27 9.34 
Sept. 8.39 8.40 8.39 8.39 8.34 8.35 
Oct. 7.71 7.70 7.75 7.95 7.90 
Dec. 6.55 7.05 6.93 


DECEMBER 1981 IR4 
TABLE 5.—Measured Consumptive Use Values 


Bermudagrass 
Bluegrass Seeded 


Wyo- Colo- 
ming” 


(3) (4) 


Collins. 
“Northglenn. 
*Tucson, Ariz. 


Measurements Arizona were taken the University Arizona Rincon 
Vista Turfgrass Research Center Tucson (8). Thirty lysimeter boxes, 
square and 0.60 deep, constructed 2.5 redwood and lined with double 
layer plastic, were used. the bottom center each box tube 
was attached and sealed allow addition removal water from the box 
desired. Each box was connected its tube large glass water bottle 
arranged Marriotte siphon. Water losses from the bottle were measured 
daily. Consumptive use from both and management regimes was 
measured. The management included more available water maintaining 
the water table consistently above that for the regime, heavier 
nitrogen fertilization, and over-seeding the warm season grasses with annual 
ryegrass for winter green grass. The ryegrass was seeded October and was 
factor until May. However, was negligible factor November through 
January. The water tables were maintained from the soil surface 
for the and from the soil surface for the regime. 
Tall fescue and St. Augustine grasses were tested only the management 
regime. Thus, they would have used still more water than indicated had 
been available. 

Measurements California were taken the University California South 
Coast Field Station Santa Ana (15). Measurements were obtained from field 
plots with water applications controlled automatically controllers activated 
tensiometers and depths. The data used herein are from plots 
having automatic irrigation and when tensiometer either reached 
cb. These were defined medium tensiometer treatments and probably 


Wyo- Colo- Cali- Cali- 
(1) (2) (5) (6) (7) 
Apr. 100 104 162 105 
May 114 124 103 203 134 
June 165 156 172 156 131 191 178 
July 183 168 152 190 169 125 177 151 
Aug. 134 124 146 153 145 127 170 124 
Sept. 113 109 136 108 141 
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Turf Grasses, millimeters per month 


Tall Fescue St. Augustine 
Tifgreen Santa Ana Zoysia 


Arizona 


not totally meet the criteria unlimited water availability, although water 
stress should have been minimal. 


The data Tables 3-5 were used calibrate the SCS (13) version the 
Blaney-Criddle formula for monthly estimates 


which the consumptive use (or E,,) for the month, millimeters; 
coefficient related temperature; crop growth stage coefficient 


for the month; and consumptive use factor. The coefficient related 
temperature 


which the mean temperature for the month, degrees centigrade. The 
consumptive use factor defined 


100 


which the monthly percentage daytime hours the year. Calibration 
the formula consisted calculating values (Table for each month using 
the averages the measured consumptive use, and temperature data over two 


= 
337 
(10) (11) (13) (15) (17) (18) (19) 
174 102 175 121 173 106 129 190 131 
204 150 197 154 191 156 105 242 175 
221 184 226 199 196 216 148 292 258 
173 141 199 154 158 174 181 248 109 217 
159 124 165 165 148 152 176 110 195 
134 116 132 127 153 140 109 164 190 
131 112 118 150 139 122 155 138 
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TABLE 6.—Crop Growth 


Bermudagrass 

Seeded 
Wyo- Wyo- Colo- Colo- 


(2) (3) (4) (5) 


Wheatland. 
Collins. 
Ariz. 


more seasons. should noted that also contains meteorological 
component and does not totally represent the crop which, this case, grass 
(6). The purpose the information presented Table (Site Descriptions) 
allow the user select values from Table from locations that most 
closely approximate their location. 

Calibrations other formulas beside the SCS (13) version the Blaney-Criddle 
were considered. The Blaney-Criddle formula, presented Criddle al. 
(4), has been calibrated using Wyoming data (8). However, the latter formula 
requires determination seasonal consumptive use factor dependent the 
length the growing season. The growing season length generally needs 
estimated, resulting small error. Many other formulas exist, but limitations 
availability required climatic data (e.g., solar radiation, wind, humidity, 
etc.) generally make these formulas impractical for use many locations. These 
formulas are value short term estimates consumptive use are desired. 
However, the present calibrations are concerned only with monthly data and 
are not meant used for detailed irrigation scheduling purposes. Doorenbos 
and Pruitt (5) have developed modification the Blaney-Criddle formula 
which permits input other variables, beside temperature and day length, without 
direct measurements. This method, however, appears have little advantage 
where data from throughout region are available for calibration the SCS 
(13) version the Blaney-Criddle formula. 


The potential evapotranspiration data from Wyoming and Colorado represent 


Month 
Jan. 1.02 0.83 1.17 
Feb. 1.03 0.71 1.07 0.44 
Mar. 1.02 0.97 1.27 0.42 
Apr. 1.72 1.76 1.16 1.28 0.83 
May 1.89 1.32 0.89 0.74 1.22 0.81 
June 1.47 1.08 1.31 1.08 0.94 0.61 0.81 0.76 
July 1.31 0.92 0.85 1.01 1.02 0.75 0.68 0.58 
Aug. 1.21 0.87 1.01 1.01 0.90 0.73 0.53 
Sept. 1.27 1.15 1.01 1.15 0.83 0.61 0.74 0.58 
Oct. 1.75 1.34 1.69 0.73 0.62 1.03 0.77 
Nov. 1.40 0.88 0.59 0.95 1.02 
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city-wide averages. Thus, the calibrated will not necessarily estimate the 
maximum potential evapotranspiration for particular plot turf. was found 
the original experiments (8) that the portion the lawn located the south, 
direct sunlight during most the daylight hours, will have 10% more E,, 
than the city-wide average. expected, that portion the lawn located 
mostly shaded areas had 10% less E,, than the city-wide average. 

The data from Arizona and California were collected from plots and are 
also good estimates city-wide evapotranspiration rates. The data from Arizona 
were collected from within city, while the data from California were collected 
experimental farm adjacent large metropolitan area. 

was possible Arizona and California grow both cool season grasses 
(Kentucky bluegrass and tall fescue) and warm season grasses (St. Augustine, 
Zoysia, Santa Ana and Tifgreen hybrid bermuda, and seeded bermudagrass). 
both locations during the warmer part the year, the cool season grasses 
had higher evapotranspiration rates than did the warm season grasses. The tall 
fescue Tucson, and Kentucky bluegrass Wheatland decreased their evapo- 
transpiration rates midsummer due heat stress. This phenomenon was 
not noted for tall fescue California. was also noted that the warm season 
grasses could mildly water stressed without loss esthetic characteristics. 
The high management regime (easier access water) was obtained maintaining 
the water table closer the surface. doubtful, however, that 


turfed area could irrigated with the precision necessary maintain mild 
water stress. 


Fescue St. Augustine 
Tifgreen Santa Ana Zoysia Cali- Cali- 
(10) (11) (13) (15) (16) (17) (18) (19) 
1.17 1.19 1.36 0.86 1.15 0.86 
0.90 0.47 1.03 0.46 1.37 0.47 1.01 0.66 0.65 0.54 
1.27 0.49 1.24 0.46 1.15 0.42 1.08 1.08 0.91 0.57 
1.38 0.81 1.39 0.96 1.37 0.84 1.44 1.50 0.86 1.04 
1.23 0.90 1.18 0.93 1.15 0.94 0.91 1.45 0.75 1.05 
0.94 0.78 0.96 0.85 0.83 0.92 1.06 1.24 0.61 1.10 
0.67 0.54 0.77 0.59 0.61 0.67 1.10 0.95 0.65 0.84 
0.68 0.53 0.71 0.62 0.71 0.63 0.95 0.76 0.66 0.84 
0.70 0.61 0.69 0.67 0.80 0.73 0.83 0.86 0.58 0.99 
0.92 0.79 0.82 1.05 0.97 0.85 0.74 1.08 0.47 0.96 
1.12 0.75 1.15 0.66 0.91 0.70 0.89 0.74 0.58 0.84 
1.02 0.45 0.75 0.65 0.74 0.60 0.86 0.79 0.73 1.07 
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Summary 


Turf grass evapotranspiration data, taken from lysimeter experiments carried 
out over several years several locations, are presented. These data have 
been used calibrate the SCS (13) version the Blaney-Criddle method for 
estimating evapotranspiration. The results should yield reasonably accurate 
estimates city-wide water requirements for turfed areas many locations 
the western U.S. The data were developed with the purpose using them 
for water supply studies and not for short-term irrigation scheduling procedures. 

warm regions, the data indicate that warm season turf grasses should 
used because they consume consideraly less water. transition regions, where 
summer soil temperatures exceed 25° users should aware that the cool 
season grasses will exhibit heat stress that reduces their evapotranspiration rates 
below what would expected with the warm temperatures. 


ACKNOWLEDGMENT 


The work upon which this publication based was supported part 
funds provided the Office Water Research and Technology, U.S. Department 
the Interior, under Public Law 88-279, the Water Resources Research Act 
1964, acting through the Wyoming Water Resources Research Institute. 


American Water Works Association, Data for Water Utilities 1970 and 
1965, American Water Works Association Statistical Report No. 20112, American 
Water Works Association, Littleton, Colo., Dec., 1973. 

Beard, B., Turfgrass, Science and Culture, Prentice-Hall, Inc., Englewood Cliffs, 
N.J., 1973. 

Cotter, J., and Croft, B., Application Practices and Landscape Attributes 
Associated with Residential Water Water Resources Research Institute 
Report No. 049, New Mexico Water Resources Research Institute, Las Cruces, 
Nov., 1974. 

Criddle, D., Harris, K., and Willardson, S., Use and Water 
Requirements for Technical Bulletin No. (Rev.), Office the State Engineer, 
State Utah, 1962. 

Doorenbos, J., and Pruitt, O., for Predicting Crop Water Require- 
Food Agricultural Organization the United Nations Irrigation and Drainage 
Paper No. 24, Food and Agriculture Organization the United Nations, Rome, 1977. 

Jensen, E., ed., Use Water and Irrigation Water 
report prepared the Technical Committee Irrigation Water Requirements 
the Irrigation Drainage Division ASCE, 1973. 

Kneebone, R., Use Proceedings Arizona Turfgrass Conference, 
Tucson, Ariz., 1979, pp. 5-8. 

Kneebone, al., Requirements for Urban Project Completion 
Report Office Water Research Technology Project B-035-WYO, Wyoming 
Water Resources Research Institute, Laramie, Wyo., Sept., 1979. 

Linaweaver, P., al., Study Residential Water Use,’’ Federal Housing 
Administration Technical Studies Program, U.S. Government Printing Office, Wash- 
ington, D.C., 1967. 

Pochop, O., al., Requirements and Application Rates for Lawns,”’ 
Water Resources Series No. 71, Wyoming Water Resources Research Institute, Laramie, 
Wyo., Feb., 1978. 

Pochop, O., and Borrelli, J., Use Kentucky Proceedings 


4 s 
4 


TURF GRASS COEFFICIENTS 


Arizona Turfgrass Conference, Tucson, Ariz., 1979, pp. 18-20. 

Quackenbush, H., and Phelan, T., Water Requirements 
Journal the Irrigation and Drainage Division, ASCE, Vol. 91, No. IR2, Apr., 
1965, pp. 11-19. 

Soil Conservation Service, Water Technical Release No. 
21, USDA-Soil Conservation Service, Engineering Division, 1967. 

Tovey, R., Spencer, S., and Muckel, C., Requirements Lawn Grass,”’ 
Transactions the American Society Agricultural Engineers, Vol. 12, No. June, 
1969, pp. 356-358. 


Youngner, B., Water Use—University California Proceed- 
ings Arizona Turfgrass Conference, Tucson, Ariz., 1979, pp. 8-17. 


IR4 341 


x 

+ 

3 


DECEMBER 1981 


PERCOLATION RATES FROM NEUTRON 
PROBE MEASUREMENTS 


Ted ASCE and Larry Mahar,’ ASCE 


INTRODUCTION 


The design irrigation systems and facilities dispose wastewater generally 
requires knowledge soil percolation rates given site. the past, soil 
percolation rates have been established with lysimeters (3,12,15,20,21,31). Three 
types lysimeters, and their characteristics, are identified Table namely: 
(1) Soil block; (2) filled-in; and (3) unconfined. Lysimeters, however, have many 
disadvantages when they are used determine percolation rates. Soil block 
and filled-in lysimeters, physically confine runoff which results unnatural 
ponding the surface. For filled-in lysimeter the disturbance associated with 
the removal and replacement the soil obviously has significant influence 
the percolation rate measured. air-water interface present, 
the case with many unconfined lysimeters, gravitational forces must overcome 
capillary tension forces before percolate can leave the soil. vacuum may 
applied through the lysimeter, but only percolate held capillary tension 
equal less than atmospheric pressure will removed. With the unconfined 
lysimeter the lack side walls precludes definition the soil block from 
which the percolate collected. Further, the soil dry, water may flow 
around the unconfined lysimeter; wet, the lysimeter may function 
drain for the surrounding soil mass. Finally, the mere presence lysimeter 
and its location and method removing percolate can cause modifications 
the natural soil moisture-tension profile and, hence, percolation characteristics. 
Several modified moisture-tension profiles which illustrate this problem are shown 
Fig. 

recognition the limitations associated with the use lysimeters measure 
quantities soil percolate and percolation rates, and part long-term 
project establish hydrologic conditions spray irrigation facility (30), 
study was undertaken investigate the feasibility using subsurface neutron 
moisture probes estimate soil percolation rates. The results from the study 
are reported herein. 

Prof. Civ. Engrg., Oregon State Univ., Corvallis, Oreg. 97331. 

Engr., ERTEC-Western, Long Beach, Calif. 

Note.—Discussion open until May 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication May 1981. 
This paper part the Journal the Irrigation and Drainage Division, Proceedings 
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TABLE 1.—Types and Characteristics Lysimeters 


Characteristics 
(2) 
Soil block Confining side walls built around essentially undisturbed block 


(monolith soil; partly open bottom installed through which percolate 
undisturbed) collected. 


Filled-in Container with vertical walls, open top, and open bottom 


through which percolate collected; filled with soil that has 
been removed from its original location and recompacted the 
in-situ condition. 

Percolate passing through in-situ undisturbed soil collected 


funnel pan placed under the soil and then conveyed the 
surface through tube. 


Unconfined 


Soil Moisture Tension 


shallow lysimeter 
with bottom drainage 


Depth (ft) 


shallow lysimeter 


evapotranspirometer sealed bottom 
with constant 


with drainage 
water table 


FIG. 1.—Soil Moisture Tension Versus Depth for Several Types Lysimeters Under 
Conditions Few Days After Precipitation (after Pelton, 1961) 
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Use Neutron Prose ror Moisture 


The development nuclear devices measure situ soil properties dates 
back 1941 when Pontecorvo (22) published paper describing device 
detect fluid bearing zones for petroleum exploration. Considerable research was 
conducted the measurement soil moisture and density nuclear methods 
during the 1950s and 1960s (2,4,5,7,8,9, which 
resulted increased confidence the use nuclear devices. Recently, 
subsurface nuclear probes have been used determine situ moisture around 
drilled shafts (6), downward movement water resulting from natural precipi- 
tation (14), the measurement situ soil moisture and density tailings 
dam (19), and changes subgrade moisture under roadbed (29). 


Connecting Cable 
Readout Device 


Preamplifier 


Detector Tube 
(with 


Collision with 
Hydrogen Atom 


Radioactive Source 
Neutron 
Bottom Cap and Shield 


(a) Probe (b) Standard (Shield) 


FIG. Neutron Moisture Probe 


The present configuration subsurface nuclear probe, specifically neutron 
probe, shown Fig The probe consists bottom cap and shield, radioactive 
source, detector tube, and preamplifier. The probe used conjunction with 
standard. The standard has dual purpose that provides the reference 
count used the calibration procedure for the probe arrive soil moisture 
content and acts shield against radiation exposure. 

obtain subsurface readings, the probe lowered into the ground through 
access (monitoring) tube. The access tube usually thin-walled aluminum 
irrigation tubing steel pipe whose diameter slightly larger than the diameter 
the probe. close fit between the access tube and the probe essential 
obtain accurate results. eliminate direct contact with the radioactive source, 
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the probe, located inside the standard, placed top the access tube 
and then lowered down the tube. 

The following events occur the determination soil moisture content with 
the neutron probe (refer Fig. 2): 


The depth moisture probe emits fast neutrons from 241 Am-Be source 
radially into the strata. 

The fast neutrons are slowed down result elastic collisions with 
the hydrogen atoms the strata. 

Some the slowed (thermal) neutrons are backscattered toward the detector 
region the probe and are absorbed the boron trifluoride (BF,). Highly 
ionized alpha particles and atoms lithium are produced upon the absorption 
the slow neutrons the 

The ionized particles and atoms lithium produce electrical pulse 
signal which amplified the preamplifier and transmitted through 
connecting cable readout device the surface. 


The detected neutrons, expressed count rate, are correlated soil moisture 
content through calibration curve. detailed description these events and 
the procedures employed develop soil moisture content calibration curve 
given Bruggers (1). 

The moisture content that determined the neutron method reported 
units density. The moisture content reflects both the amount free and 
held water the soil mass. (Free water that water which may move through 


soil mass under the influence gravity while held water that water which 
does not move under the influence gravity.) 


Prose Measurements 


Consider quantity water applied the surface soil column which 
subsequently percolates through the column shown Fig. 3(a). the water 
percolates past representative point the soil column two timedependent 
relationships may observed. shown Fig. the moisture content, 
starting from initial value, will increase (leading edge) then decrease (trailing 
edge) with increasing time. there increase storage the soil, the 
moisture content will return its initial value. shown Fig. 3(c) the velocity, 
zero, also increases and decreases with increasing time, eventually 
returning zero. Referring Fig. noted that any time the 
amount water moving through the soil column may assumed equal 
the change moisture content the soil column, that is, Amc,. The amount 
water the soil column moving velocity noted Fig. 3(c). 
inspection, the percolation rate time PR,, may calculated from 
PR, 


which, y,, density water. 
determine percolation rates employing Eg. is, therefore, necessary 
determine the moisture content and velocity relationship with time given 
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point soil column. practice, possible determine the moisture 
content versus time relationship given point with subsurface neutron 
probe, but not possible determine the velocity versus time relationship 
given point. However, the velocity versus time relationship can determined 
over representative volume material follows: 


Applied 
Woter 


| 


Moisture Content 
Velocity 


| 4 ‘ 
bitrory 
or 
point - 


oO 


20 25 3% 
Elapsed Time (minutes) 


FIG. 4.—Soil Moisture Content Versus Elapsed Time Two Levels 


Determine the soil moisture content versus time relationship two levels, 
and II, soil column separated distance (refer Fig. 4). 

Divide the ordinate values moisture content the respective peak 
moisture contents obtained the two levels; the result two normalized moisture 
content versus time relationships approximately the same shape but with 
peak values different times. 


Superimpose the two normalized moisture content versus time relationships 


| 
(a) 
FIG. 3.—One Dimensional Percolation Soil Column 
7 
evel I 
Jevel 
5 
4 
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Calculate the velocity the water function normalized moisture 
content noting the elapsed time between equal normalized moisture contents 
the leading and trailing edge, respectively, level and and dividing 
the distance between the levels, the elapsed time; the result the relationship 
shown Fig. 

For either level time determine the normalized moisture 
content from Fig. and the corresponding velocity, from Fig. for the normalized 
moisture content time associated with the leading trailing edge, respectively. 
considering several values time the velocity versus time relationship may 


Normolized Moisture Content 


Elapsed Time (minutes) 


FIG. Normalized Moisture Contant Versus Elapsed Time Two 
Levels 


— 


leading edge 


trailing edge 


Normalized Moisture Content 
FIG. 6.—Velocity Versus Normalized Moisture Content 


The validity the approach outlined above based primarily the assumption 
that the velocity the water can determined from knowledge the 
elapsed time between equal normalized moisture contents. the authors’ 
knowledge there precedent for this assumption for unsaturated flow. For 
saturated flow with changes soil fabric during flow, Darcy’s Law provides 
validity for the procedure. Specifically, Darcy’s Law may written 


wailing edge 
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Elapsed Time (minutes) 


FIG. 7.—Velocity Versus Elapsed Time 
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Soil Moisture Content (pcf) 


Elapsed Time (minutes) 


FIG. Soil Moisture Content Versus Elapsed Time 
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The assumption stated would mean the calculated velocity simply associated 
with the average permeability and average hydraulic gradient between two levels. 
many applications average percolation rate required rather than the 


Moisture Content 
Moisture Content 


time time 
(a) Peaking (b) Leveling 


Shapes for Soil Moisture Content Versus Time 


Access Tube 


Wire Screen 
covered with 
Filter Paper 


FIG. Test Setup for One-Dimensional Percolation Study 


percolation rate specific time. The average percolation rate may determined 
follows: 


Divide the moisture content versus time curve level into increments 
time shown Fig. the increments should smaller the nonlinear 
regions; they can larger the linear regions. 


Inflow 
Gravel 
Level 
— 
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For each time increment estimate the average change moisture content 
and the average normalized moisture content. 

Estimate the average flow velocity each element with knowledge 
the velocity versus normalized moisture content relationship determined for 
the test (refer Fig. 6). 

Calculate the percolation rate for each time increment using equation (1). 

Calculate the amount water (in units length) associated with each 
time increment multiplying the percolation rate for each time increment 
the increment time. 

Sum the amounts water (in units length) associated with each time 


increment and divide the total amount water the total elapsed time for 
the event. 


The soil moisture content versus time relationship may represented 
either the two characteristic shapes shown Fig. small quantity 
water applied over short time interval relationship between 
soil moisture content and time results [Ref. Fig. 9(a)]. large quantity 
water applied over short time interval and the peak moisture content equals 
the capacity the soil, small quantity water applied over long 
time interval, relationship between soil moisture content and time 
results [Ref. Fig. 9(b)]. The procedure presented applies either 
relationship between soil moisture content and time. 

Finally, interest note that the procedure presented determine 
soil percolation rates based differences moisture content. Therefore, 
systematic errors associated with the determination the absolute value 
moisture content with the neutron probe will cancel out and not affect the 
determination the percolation rate. 


the laboratory program evaluate soil percolation rates column 
well-graded sand was used represent one-dimensional model. The test setup, 
shown Fig. 10, consisted two gal (0.21 drums containing 4.1 
(1.3 column sand. aluminum access tube was centered the drums 
and the sand was placed around the access tube pouring through #30 
sieve maintained constant drop height (300 mm) over the surface. 
in. (100 mm) layer gravel was placed the base the column. The 
gravel allowed uniform drainage over the cross section the column. allow 
water applied rate less than the flow capacity the sand, wire 
screen covered with filter paper was placed top the column. Finally, 
thin layer gravel was placed top the filter paper and screen 


minimize disturbance the sand during application the water. The procedure 
conduct test was follows: 


Initial moisture content readings were taken three levels: 1.1 ft, 2.1 
ft, and 3.1 (0.35 0.66 and 0.97 m), shown Fig. 10. 

With the subsurface neutron probe level water was applied constant 
head until the total amount water associated with the test had been added. 
The water level was then allowed fall. 


Change in 
Moisture 
Content, in 
pounds per 
cubic foot 
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TABLE 2.—Example Calculation Percolation Rates 


Elapsed Time 
between Equal 
Normalized Percolation 
Moisture Velocity, in Time at Rate, in 
Contents, inches per Level |, inches per 
in minutes minute in minutes minute 


Normal- 
ized 
change 
in mois- 
ture 
content 


0.100 


0.090 


0.080 


0.070 


0.060 


0.050 


0.040 


0.030 


Percolation 


0.010 


0.0 


Elapsed Time (minutes) 


FIG. 11.—Percolation Rate Versus Elapsed Time 
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: Level ing ing ing ing ing ing ing ing 
4 ! edge edge edge edge edge edge edge edge 
(1) (4) (5) (6) (7) (8) (9) (40) (11) 
0.52 0.37 0.1 4.21 240 2.85 0.05 2.6 os 0.024 0.0004 
1.04 0.74 0.2 5.45 109 2.20 0.11 3.6 _ 0.037 0.002 
’ 1.56 1.1 0.3 6.42 7S 1.87 0.16 4.2 62 0.047 0.004 
z 2.08 1.48 0.4 7.10 54.5 1.69 0.22 4.6 34.5 0.056 0.007 
2.60 1.85 0.5 7.74 37.5 1.55 0.32 5 24.5 0.065 0.013 
‘ 3.12 2.22 0.6 8.22 28.6 1.46 0.42 5.3 20.0 0.073 0.021 
3.64 2.59 0.7 8.57 20.0 1.40 0.60 5.6 15.6 0.082 0.035 
4.16 2.96 0.8 8.82 16.0 1.36 0.75 6.0 13.8 0.091 0.050 
4.68 3.33 0.9 9.60 13.48 1.25 0.89 6.5 10.5 0.094 0.067 
$.20 3.70 1.0 11.43 11.43 1.05 1.05 7.4 74 0.088 0.088 
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timer was started the instant the water was applied and moisture 
content readings were taken every sec level When increase moisture 
content was observed the probe was lowered level and readings were 
taken sec intervals. When increase was observed level the probe 
was moved back level where another sec reading was taken. Thereafter 
the probe was moved consecutively all three levels sec intervals and 
readings were taken. 

Readings were taken until either the count ratio was equal the initial 
reading until the moisture content versus time curve became linear. 

The percolate was collected the base the column and weighed. 


The data obtained from the test can represented plot moisture 
content versus time for each the three levels. The data was analyzed 
estimate the percolation rate versus time relationship using the procedures 
discussed the preceding section. example calculation given Table 
for the event shown Fig. recorded for levels and Figs. 
and support the example calculation presented. summary the test results, 
the percolation rate versus time relationship, shown Fig. 11. 


Test 


The validity the proposed procedure determine percolation rates was 
assessed calculating the total amount percolate for given test and comparing 


TABLE 3.—Example Calculation Amount Water Percolating Through Soil Column 


Average 


change 
moisture Average Perco- 
content, normal- lation 
Increment pounds ized Velocity, Amount 
time per cubic moisture inches per inches water, 


(2) 


anne 


WN 


Un oo 


in. 25.4 mm. 


Ele- 
ment 
num- 
ber 
3.0 0.29 1.9 
2.9 1.5 2.3 
4.7 1.25 1.5 
1.00 1.00 1.4 
3.0 .017 2.2 
1.0 .002 3.1 
otal 
22.7 
22.0 
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the amount collected. The calculation the amount percolate requires 
the use the percolation rates. Therefore, favorable comparison between 
the amount percolate calculated that collected would substantiate the 
intermediate calculation the percolation rate. The total amount percolate 
may determined the procedure given preceding section estimate 
the average percolation rate. deviation the procedure required the 
last step follows: 


Sum the amounts water (in units length) associated with each time 


increment and multiply the cross-sectional area the drum and the density 
water. 


Best fit line, 


— 


Weight Water Collected 


FIG. 12.—Weight Water Versus Weight Water Collected Level 


example calculation given Table for the test results presented 
Fig. Fig. supports the example calculation presented. 

total laboratory tests were conducted with the amount water applied 
ranging from (4.5 kg). The percolation rates and total amount 
percolate were determined for each test using the data from levels and 
II, and and and III. Fig. shows the comparison between the total 
amount percolate calculated using the data from levels and III and total 
amount collected. (Note: Level measurements are missing for one the 
tests, thus, only nine data points are shown.) The least squares best fit line 
for the data set, with correlation coefficient 0.98 also shown. The slope 
the regression line close the 1:1 line and suggests that the estimated 
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percolation rates are very reasonable compared the measured values. 

Fig. shows the comparison between the total amount percolate calculated 
using the data from levels and and and III. The least squares best fit 
lines for the data sets, with correlation coefficients 0.96 (level and II) 
and 0.98 (level and III) are also shown. The comparisons are not favorable 
that shown Fig. 12. The relationship obtained Fig. may explained 
considering the sphere influence the neutron probe. van Bavel, 


al. (28) determined approximate relationship for the radius, the sphere 
influence follows: 


355 


100 
moisture content percent volume 


For the tests conducted the radius the sphere influence between 
in. (460 mm) low moisture contents and in. (230 mm) high moisture 
contents. Therefore, the low moisture contents associated with the start 
test, for the level readings, the radius influence would extend above 
the top the soil column. When the water reservoir was created above the 
soil column conduct test, the zone influence would extend into the 
reservoir and result higher moisture content readings than exist the soil. 
Calculated percolation rates based level data would, therefore, higher 
than actually occurred. Also calculated weights percolate water would 
greater than the amounts percolate water collected shown Fig. 13. 


The procedure estimate soil percolation rates with subsurface neutron probes 
was presented preceding section. The procedure applies the determination 
percolation rates for both laboratory and field conditions. The procedure 
employed determine soil percolation rates situ will differ from the procedure 
employed the laboratory two aspects, namely: (1) The installation the 
access tube; and (2) the time increments which moisture content readings 
are taken. 


Installation the access tube should satisfy three requirements (29): 


There should snug fit between the access tube and surrounding soil. 
The access tube should watertight prevent water and moisture from 
entering. 


The access tube should straight and circular allow passage the 
probe. 


least three methods are commonly employed meet these requirements 
follows: 


Driving the access tube closed-end means drop hammer (applicable 
only very soft soil deposits). 


Driving open-end access tube and removing the material from inside 
the tube with auger. 
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Pre-augering hold and driving the access tube closed-end into the hole. 


would appear that any one these three techniques would suitable. 
The time increment which moisture content readings are taken will depend 
great degree the soil conditions given site. Fig. shows the 


Best fit line, 
Level 


\Best fit line 
Level 


@ 


200 240 280 320 
Weight Water Collected 


FIG. 13.—Weight Water Calculated Versus Weight Water Collected Levels 
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moisture content versus time relationship obtained level (0.6 
following the spray application spray irrigation facility. The soils the 
site consist silty fine medium glacial outwash sand overlain in. 
(80 mm) layer organic topsoil and vegetal cover. The time scale for the 
moisture content relationship for the materials the site hours, owing 
the relatively low permeability the soils the site. Moisture content 
readings were taken min intervals. Further, considerable time period 
was required obtain the complete moisture content relationship this particular 
site. 

The results from the laboratory test program suggest that moisture content 
readings should not taken close the surface the soil deposit and the 
distance between the levels which readings are taken should (0.31 
less. Additional research work needed more accurately define 
appropriate distances from the surface and between successive levels function 
soil type and moisture content. 

van Bavel, al. (28) suggest that the neutron probe can lowered into 
access tube and the density detected neutrons, expressed count 
rate, can observed function depth. The count rate will increase 
with depth until the soil-air boundary outside the sphere influence 
the probe. this depth the count rate will level off and the readings are 
associated with body. This concept could easily applied 
field situation determine the minimium reading depth. 

Finally, the conclusions drawn the paper are associated with moist 
homogeneous column well-graded sand. the field layered soil profiles 
(in which the materials have significantly different permeabilities) shallow 
ground-water table conditions can exist. The validity the proposed procedure 
under these conditions questionable and requires further evaluation. 


Summary 


procedure presented estimate soil percolation rates with subsurface 
neutron probes. The procedure based upon the soil moisture content versus 
time relationship obtained two levels soil column soil profile which 
one-dimensional percolation occurring. The results from laboratory tests 
employing column well-graded sand which percolation was occurring 
were used validate the procedure. The calculated quantities water for 
the percolation tests compared favorably the actual amount water 
collected. Therefore, the procedure proposed appears acceptable for 
homogeneous deposits well-graded sand. The validity the procedure for 
layered soil profiles shallow ground-water table conditions must investigated. 
The procedure proposed simpler than traditional methods (e.g., lysimetric) 
measure percolation rates. 

When using subsurface neutron probes estimate percolation rates situ, 
consideration must given the method installation the access tube, 
the time intervals which moisture contents will determined particular 
site, the depth from the ground surface the level where moisture content 
readings are taken, and the distance between levels where moisture content 
readings are taken. 
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DIMENSIONLESS STREAM ADVANCE 
SLOPING BORDERS 


INTRODUCTION 


mathematical model (7) border irrigation has been developed and found 
both reliable and sufficiently inexpensive execution comprise practical 
design tool. the other hand, does require from fraction minute 
several minutes execution time digital computer intermediate speed 
and memory capacity (e.g., Burroughs 6700) simulate the advance, recession, 
runoff, and ultimate distribution infiltrated water stemming from given inflow 
and application time and bed slope, roughness, and infiltration. Furthermore, 
the results are course valid only for that one set input parameters. 
common with results numerical methods general, the solution obtained 
gives clue regarding solution behavior general. Behavioral trends can 
noted only through repeated solution using many different combinations input 
variables. Because the large number variables governing the outcome 
any one irrigation, comprehensive examination the effect each, 
dimensional terms, nearly impossible. Presentation the results alone would 
require pages and pages graphs. 

The problem discerning the general behavior solutions greatly simplified 
the use dimensionless notation. Use such notation reduces the number 
independent parameters governing the solution tractable value. The result, 
several pages dimensionless graphs, displays the overall behavior the advance 
function under the most varied input conditions found practice. Furthermore, 
through interpolation between curves, the advance for any practical combination 
input variables can obtained directly, without the use electronic computing 
equipment. 

While the complete mathematical model outlined (7) has the capability 
predicting the entire irrigation, including recession and runoff, the dimension- 
less solution graphs herein are restricted the advance phase alone. Recession 
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controlled more parameters than advance, and representation the complete 
theoretical spectrum recession curves not yet possible, even with full 
use dimensionless notation. This problem currently under study, and the 
writers express their optimism, regarding satisfactory resolution the future. 

initial study dimensionless advance was presented Ref. pointed 
out Ref. the range variation the dimensionless variables employed 
excluded some cases practical importance. Furthermore, the choice 
dimensionless parameters was not optimal and required excessive number 
cycles presenting the data logarithmic scales. Finally, the response 
the advance function cutoff the inflow was ignored. All these questions 
are treated the present paper. 


The governing equations imply conservaton mass and equilibrium forces 
elements the surface stream. The forces can assumed balance, because 
accelerations are negligible border-irrigation flows (3,6). The upstream boundary 
condition the discharge hydrograph, typically, inflow constant, until 
cutoff the end the application time, The downstream boundary condition 
during stream advance zero depth and discharge the wave front. Initially, 
depth and discharge are zero everywhere the border. 

The mathematical expressions these conditions the following: 


2 2 


x=x,<L 


Eqs. 1-5, discharge per unit width; unit discharge crossing the 
face element the surface stream; constant unit inflow; 
surface-flow depth; subsurface depth (volume infiltrated per unit plan area 
border); distance down the border; length the advancing stream 
(assumed for this study always less than field length, L); time; and 
cutoff application time. The subscript, refers the upstream end 
the border; refers the upstream face element the surface stream 
and the downstream face. The length the element while 
increment time. The bar over variable signifies time average during 
the tilde represents distance average over The bottom slope given 
S,; represents the resistance (drag) the bed and vegetation per unit 
plan area surface stream. related the so-called friction slope, 
d/y, and generally function, y), depth and discharge section, 
for any given soil-surface and crop geometry. The depth infiltration any 
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(x), that location; (x) the inverse the advance function, 
(t), and represents the time the stream advanced station 


Equations 1-5 can put dimensionless form first introducing nonzero 
reference variables and For the moment, these are undefined, 
but possess time, length, length, and length, respec- 
tively. Then dimensionless variables are defined follows 


and are substituted into the governing equations. When divided 
XY, and Eq. Y’, the following dimensionless forms result 


inwhich 


* 


The quantity, ratio local bed slope, bed slope, the head 
end the border; plane border The term, ratio local 
drag per unit plan area, drag per unit area, computed with reference 
depth and discharge, y)/D; and d(Q, Y). The infiltration depth 


it=—; x*=—; z*=—; 
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dimensionless time and appear Eq. instead dimensional time, 
and the reference infiltration depth, must defined function 
the characteristic time, Z(T), such that 


The functional form usually, but need not be, the same the functional 
form The expression terms depends upon the exact form 
the actual and reference infiltration functions; this question dealt with 
later paragraph. 

The solution Eqs. and plane border controlled the dimensionless 
parameters, ¥*, K*, and D*, the dimensionless boundary values, 
and the form the dimensionless infiltration and resistance functions, 
and d*. 

Selection Dimensionless Governing Parameters.—Substantial simplification 
these equations achieved judicious choice the reference variables 
that and was set normal depth for the given inflowing discharge, 
slope, and roughness. This led the results, S$. The reference infiltration 
depth was defined terms the reference time, which was 
the given infiltration function. Ref. was determined the condition 
Then was found from the additional desired condition 
Furthermore, Ref. the two-parameter Kostiakov infiltration formula, 
which and are constants, was used. result, the system 
Eqs. 7-11, with and indefinitely large, was governed only two 
parameters, and and dimensionless-advance solutions were presented 
two-parameter family graphs. Clemmens (2) pointed out that for given 
range the physical variables, and Manning the same information 
could presented over smaller range x*, requiring instead 
and the result was again two-parameter family graphs and 
were still considered larger, respectively, than any computed and x*), with 
the governing parameters. 

Ref. the two-parameter family could have been reduced one governed 
single parameter, because level borders, with only, were considered. 
With and ¥*, K*, and D*, each set unity define and 
and the two remaining solution-governing parameters would 
the sole governing parameter. the other hand, the requirement 
for all soils, however tight permeable they may (e.g., for the entire range 
Kostiakov and found practice), leads enormous variation 
and, thus, From this stems the need for many, many cycles logarithmic 
paper which present graphical advance data. Furthermore, because advance 
(and recession) after cutoff also was considered Ref. proved 
solution was governed and 

The question arises: how may the dimensionless governing parameters 
¥*, etc., can controlled independently the physical 
parameters the irrigation? corollary question is—what the minimum 
number parameters governing the solution Eqs. 7-11? 
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Eqs. 12-15 and the last two parts Eq. can rewritten equations 
defining the interrelationship between the reference variables and the dimension- 
less governing parameters terms the physical parameters the irrigation 


The six equations, Eqs. 17-22, must satisfied, whatever are the values ultimately 
taken the six dimensionless parameters, ¥*, K*, D*, and 
the four reference variables, and other criteria need met. 
Thus, there considerable leeway assigning values these quantities. 
Comparison the number variables given values, namely 10, with 
the number equations relating these values, six, suggests that four variables 
can chosen arbitrarily, with the remaining six stemming from solution 
the six equations. the other hand, the particular structure these equations 
introduces constraints and reduces the degree freedom that can exercised 
choosing those variables assigned arbitrary values. For example, Eq. 
19, alone, shows that even the three quantities, and cannot given 
arbitrary values. 

Thus, mentioned previously, Ref. was set normal depth for 
and S,, and ¥*, and were each set unity. This left 
and govern the dimensionless solution together with whatever parameters 
were introduced the infiltration and friction-slope functions. Ref. 
was assumed larger than the largest considered, that the last analysis, 
(like L*) played role. 

were set unity. The dimensionless solution was then governed and 
and parameters from the infiltration and friction functions. Again, 
was set excess all reduce the total number governing parameters. 

set unity, and Eqs. and yielded and Equation gave 
which along with the parameters the infiltration and friction functions, governed 
the dimensionless solution. 

Infiltration Functions.—In Refs. and the infiltration function adopted 
was the Kostiakov formula 


Z(T) 
Xd(Q, 
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The same functional form was used for This led simple result for 


Thus, addition whatever other dimensionless parameters governed the 
solution, the parameter, stemming from the infiltration function also played 
role. 

the variable representing irrigation times becomes quite large, the Kostiakov 
formula gives too small infiltration rate, and fails reproduce the constant 


rate exhibited many soils large infiltration times. The following modification 
the Kostiakov equation 


allows nonzero final infiltration rate, but now additional governing parameter 
arises, can seen from the corresponding dimensionless infiltration function 


obtain Eq. 26, the functional forms and were again made identical, 


The two parameters stemming from the infiltration formula 
are now and b*, the latter given 


The Philip formula (5), special case Eq. 25, and 
leads but one infiltration-based free parameter, given 


possible develop more direct approach than Eq. for synthesizing 
infiltration function which reflects both the high initial rate characteristic 
the Kostiakov formula and constant final infiltration rate two-branch 
function suggested, the branches meeting the time, t,, when the rate given 
the Kostiakov formula equals i,. Then 


a 


Note that 


Such function more easily obtained from plot logarithmic paper 
measured infiltration rate against time than the function Eq. 25. were 
again chosen equal dimensionless infiltration, z*, function 
dimensionless infiltration time, would have four branches depending upon 
function can described with but two branches, follows, the reference 
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1/(a—1) 


The two-branch function, Eq. 31, then, introduces two additional dimensionless 
parameters, and into the set governing solution Eqs. 7-11. 

Resistance Functions.—Next, solution-governing parameters stemming from 
the choice resistance function are considered. With the Manning formula, 
the drag per unit area given 


SI, system units; 1.4859 ... the English system), and 
the Manning roughness coefficient. Now, the reference drag per unit area 


Thus, additional parameters are introduced into 7-11. 
With the Sayre-Albertson formula (4), the drag per unit area 
2 


2 


investigated irrigation borders Kruse (4). With d(Q, Y), the relative 
drag per unit area given 


. 2 
log 
Evidently, this form the resistance formula introduces additional solution- 
governing dimensionless parameter x*, defined 
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Curves 


From the foregoing analysis, clear that the minimum number dimension- 


w 
a 


OIMENSIONLESS DISTANCE 


FIG. 1.—Dimensionless Advance; 0.0 


DIMENSIONLESS TIME 


OIMENSIONLESS DISTANCE 


FIG. 2.—Dimensionless Advance; 0.1 


IR4 DIMENSIONLESS STREAM ADVANCE 369 


less-solution parameters for advance smaller than field length and for advance 
times smaller than application time, one, plus those stemming from the particular 
choice infiltration and resistance formulas. With the Kostiakov and Manning 


DIMENSIONLESS TIME 


OIMENSIONLESS DISTANCE 


FIG. Advance; 0.2 


w 
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n 
z 
a 


DIMENSIONLESS DISTANCE 


FIG. Advance; 0.3 


equations, one additional parameter, introduced. Kostiakov formulation 
modified yield constant final infiltration rate per Eqs. requires 
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one additional free parameter while the Sayre-Albertson resistance formula 
introduces one more. Introduction cutoff time pertinent variable leads 
additional dimensionless solution parameter. 


100. 


DIMENSIONLESS TIME 


OIMENSTONLESS DISTANCE 


FIG. Advance; 0.4 


z 
= 


DIMENSIONLESS DISTANCE 


FIG. Advance; 0.5 


Generalized advance curves are presented for sloping borders indefinite 
length with resistance expressible terms the Manning and infiltration 
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characteristics defined the Kostiakov power law. Because finite application 
times are considered, the minimum number solution parameters three. 
principle, one-parameter family dimensionless advance curves can 


100. 


DIMENSIONLESS TIME 


DIMENSIONLESS DISTANCE 


FIG. 7.—Dimensionless Advance; 0.6 
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DIMENSIONLESS DISTANCE 


FIG. Advance; 0.7 


shown one sheet graph paper; each curve, would correspond 
different value the parameter. two-parameter family requires booklet 
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pages (as Refs. and 3), each page corresponding different value 
the second parameter. three-parameter family requires then, general, 
series booklets, each identified the value the third parameter. 


DIMENSIONLESS TIME 


DIMFMSIONLESS DISTANCE 


FIG. 9.—Dimensionless Advance; 0.8 


DIMENSIONLESS TIME 


DIMENSIONLESS DISTANCE 


FIG. 10.—Dimensionless Advance; 0.9 


Such would have been the requirement for presenting generalized advance 
data, including advance after cutoff, the otherwise desirable choice was made 
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reasonable upper limits) would cover all possible physical variation (Kostiakov) 
infiltration and bed slope, including impervious soil and horizontal borders 
special cases. But order limit presentation data essentially one 
the possibility including the series disallowed, and 
advance data for given relatively crude form, without the detail 
evident the advance curves for The latter information given 
separate set curves. 

and along with Manning and Kostiakov the respective 


w 
o 
2 
= 
Zz 
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DIMENSIONLESS DISTANCE 


FIG. 11.—Dimensionless Advance; 1.0 


resistance and infiltration functions, the characteristic depth, distance, and time, 
respectively, are 


@.1 4. 10. 100. § 
qin 
n 
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The dimensionless solution-governing parameters are then K*, The 
graphs are organized around these three parameters follows. Dimensionless 
advance, prior cutoff presented function dimensionless time, 
t*, with dimensionless infiltration, K*, and exponent, parameters. This 
effectively eliminates, for the present, governing parameter because, 
The solution curves were obtained means the numerical methods described 
Ref. 

The range dimensionless variables presented stems from the extensive 
survey dimensional irrigation parameters encountered practice described 


o 
Zz 


DIMENSIONLESS DISTANCE 


FIG. 12.—Dimensionless Advance over Impervious Bed 
Ref. The ranges adopted are 
L/sm 


100< 500m 


These lead, with the aforementioned system for nondimensionalization, the 
following ranges dimensionless time and distance 


7 
1000. 
100. 
0.10 10. 100. 1000. 
(45) 
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and dimensionless infiltration 


for a=0.0 
for 
for a=1.0 


Considering: (1) The values and resulting from input-variable 
combinations actually found the field (as opposed artificial combinations 


(Solid curves) 10. 


FIG. 13.—Maximum Advance Front-End Recession) after Cutoff; Distance 
and Time; 0.1 


extremes taken from different sets field data); and (2) the aforementioned 
limits and found from such artificial combinations, the following 
limits are expected cover the range conditions pertinent sloping borders: 


100 


100 
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The data displayed Figs. are limited these ranges. 

Advance Impervious Bed.—Examination these graphs for 
shows, not surprisingly, that this curve independent the Kostiakov 
Furthermore, the bow this curve, concave downward, has theoretical 
explanation. the early stages advance, when very small, the depth 
gradient the surface stream very large everywhere. Eq. with 


(dashed curves) 


FIG. 14.—Maximum Advance (Incipient Front-End Recession) after Cutoff; Distance 
and Time; 0.3 


clear that the middle term negligible when very small. The solution 
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Eqs. 9-11, and then the same one obtained from 7-11 
with This solution shown Fig. 12, and indeed seen 
the limit the general solution when very small. proven Ref. 
this curve has the constant slope 16/13 when plotted logarithmic paper. 
the other extreme, when very large, the advance rate becomes uniform 
the profile assumes its ultimate shape. Then, the slope the advance curve 
logarithmic paper simply unity, thus, explaining the observed curvature 


(dashed 


10. 


FIG. 15.—Maximum Advance Front-End Recession) after Cutoff; Distance 
and Time; 0.5 


the general advance curve with Further, with the particular system 
nondimensionalization employed, the dimensionless speed advance also 
unity. Indeed, under ultimate conditions, with fixed profile shape moving 
down the border, advance rate equals the velocity flow the 
Advance After Cutoff.—Once the irrigation stream has been cutoff, the 
slows. With sloping impervious bed, advance continues indefinitely 
though ever-decreasing rate. With infiltration, however, the advancing 
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Stream eventually halts and then recedes back the border. The curve 
advance after cutoff approaches some maximum advance distance, 
time, which represents the termination advance and the beginning 
front-end recession. The maximum advance distance and associated time 
are shown Figs. 13-17 for 0.1-0.9 steps 0.2. Both time (dashed 
lines) and distance (solid lines) are shown relative the values each 
cutoff. For high values K*, front-end recession occurs almost immediately 


after cutoff and the stream advances only short distance beyond 


FIG. Advance (incipient Front-End Recession) after Cutoff; Distance 
and Time; 0.7 


usually undesirable for the field length greater than the maximum 
advance distance since part the border would receive water all. But 
knowledge the maximum advance distance and its associated time can 
used construct approximate curve advance after cutoff, even the 
theoretical maximum advance greater than field length. After calculating the 
appropriate values K*, and L*, Figs. 1-11 are used with 
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the appropriate value construct the curve for advance prior cutoff. 
Then Figs. 13-17 are used determine the maximum advance distance and 
its associated time. The curve for advance after cutoff can then sketched 
through the two known points and with the 
requirement that tangent the curve for advance prior cutoff the 
time cutoff and tangent the vertical const) the end advance. 
any event, the advance curve constructed beyond simply ignored. 


(dashed curves) 


FIG. Advance Front-End Recession) after Cutoff; Distance 
and Time; 0.9 


Example 

Determine the advance curve for the following situation: 1200 (366 

From Eq. (with units meters and seconds) 
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2430 sec 40.5 min 0.675 


Gin * 7.432 ¢/s-m 

L *366m Max. advance 
= 5.08cem/hr® distance and 

a *07 time from Fig /6 

n 

Sq * 0.00) 

8O min 


Approximated 


advance curve 


Field length —~— 


Advance curve 
from Fig.8 


100 200 
DISTANCE (m) 


FIG. Advance Dimensional Terms 


250 The advance curve cutoff can choosing 
values and the corresponding The maximum possible advance 
corresponding time 2.5 min 200 min. The entire 
advance curve can now constructed Fig. 18; the curve after cutoff 
tangent advance cutoff and the vertical the maximum advance 
distance. The total advance time for the given field length from Fig. 
roughly 136 137 min, which very close the time computed with the 
computer model directly, namely 136.1 min. From the computer solution 
this example (which includes recession), the application efficiency (net applica- 
application) for in. (7.62 cm) net application was 78% with 
runoff. 


0.1344 
The advance curve can constructed from Figs. and which 
200 
160 
120 
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important note that interpolation between curves can subject 
error. However, errors estimating maximum advance distances and times 
will produce minor errors unless the field length very close maximum 
advance distance. The major source error will normally determining 
advance prior cutoff. interpolation questionable, plot advance distances 
any given time, for the values presented, and graphically determine 
the advance distance for the desired value. 


Summary 


The simplest case general border-irrigation advance characterized 
plane border length greater than any advance distance interest, and 
constant inflowing discharge prior cutoff. this case, the advance 
curve controlled the border slope, roughness, and infiltration, and 
the inflow discharge and cutoff time. the Manning formula used describe 
the roughness and the Kostiakov formula the infiltration, the advance function 
significantly reduced—to three—when the problem viewed dimensionless 
terms. advance prior inflow cutoff only interest, two dimensionless 
parameters govern the dimensionless advance curves. 

and are chosen those governing parameters (by setting 
needed cover all cases practical interest. The case impervious bed 
cannot represented all. the other hand, possible show directly 
from the curves Ref. how kinematicwave conditions 
are approached large This information useful establishing the 
range application kinematic-wave theory computation border-irrigation 
flow, procedure simpler than the zero-inertia (equilibrium) analysis employed 
for the general case. 

disallows simple comparison between zero-inertia and kinematic-wave 
solutions, but permits the whole spectrum dimensional parameters encountered 
practice represented data covering just few cycles logarithmic 
paper. 

With pertinent controlling variable, the number dimensionless 
governing parameters increases three. the cutoff time were made the 
characteristic time, dimensionless solution would governed and 
This would have the advantage allowing representation all slopes, 
including zero slope, and all infiltration rates, including zero, the one set 
curves. This entails, however, major disadvantage, the difficulty 

this work, compromise was achieved using normal depth the 
characteristic depth and setting This longer permits consideration 
horizontal borders, but these were extensively studied Ref. Prior 
cutoff, the advance curves are controlled and K*, and this two-parameter 
family presented Figs. 1-11. Advance after cutoff viewed approximately 
Figs. 13-20, which give the end points this portion the advance curve 
and some indication how the curve may sketched between the starting 
point, and the end point, 


7 
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Taken together, these curves represent all practical cases advance 
irrigation stream over plane border with inflow that constant prior cutoff. 
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RELATIVE RUNOFF CURVE 
NUMBER NOMOGRAPH 


INTRODUCTION 


many drainage analyses, the Soil Conservation Service (SCS) curve number 
(CN) runoff volume estimate may serve initial basis for design engineering 
structures (6,7). While there are often variety runoff-estimation techniques 
available having sounder basis theory than the procedure, these methods 
are more often appropriate for subsequent, data-intensive studies (2,9). 

The method weighs runoff-influencing effects watershed condition and 
has been pragmatically incorporated into variety engineering procedures 
(1,3,4,5,8). The SCS has cataloged soils hydrologic families and reported 
the variation with antecedent moisture, land condition and use, vegetative 
type and cover, and soil type. SCS soil survey studies often comprise the majority 
compiled data available for preliminary hydrologic evaluations (7,10). Whereas 
other methods may more accurately predict runoff hydrograph peaks and times 
concentration, the method can rapidly and efficiently appraise the 
volumetric sensitivity the hydrograph changes watershed condition. 

Relative increase runoff paramount technical issue urbanizing areas 
where land use changes are only allowed increased runoff can retained. 
suggested that while more theoretical statistical methods should generally 
used estimate the time-dependent magnitude discharge, analysis 


may used explore the relative dependency this magnitude watershed 
condition. 


Discharge volume may estimated the SCS dimensionless standard 
formulation (6) 


P+0.8S 


which runoff inches (cm); precipitation inches (cm); and 
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potential maximum retention inches (cm). the inch system, the parameter, 
estimated 


Relative runoff can defined the ratio volumetric runoff from watershed 
alternative hydrologic conditions. Hydrologic state indicated the effective 


20 10 ‘2 | 05 02 
. 


= 


Q,/0, = 1.00 


| 


FIG. 1.—Relative Runoff Nomograph 


CN. The definition can broadened the assumptions associated with unit 
hydrographs are employed. These standard simplifications indicate that for 
given-duration storm time concentration and runoff duration are independent 
the hydrograph ordinate. While discharge varies with excess precipitation, 
the relative shape the hydrograph does not. Thus, relative runoff can also 
taken approximate ratio peak flows. The same caution should 
taken with this latter employment would taken synthetic unit hydrograph 
derivation between storms significantly different magnitudes. 

Employing the subscript denote watershed state and runoff before 
hydrologic modification and the subscript after hydrologic modification, Eqs. 
and can combined: 


(2) 
— CN . . . . . . . . . . . . . . . . . . . 
which empirical curve number. the metric system 
2540 
CN, 
a 398 + + 
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which Q,/Q, represents relative runoff. 

should noted that unlike relative discharge estimates the rational 
formula, this ratio depends precipitation also. Abstractions from light rainfall 
influence runoff larger degree than abstractions from heavy rainfall. 

Eq. may graphically solved Fig. Vertical and horizontal nomographs 
are used intersect the left-hand and upper sides, respectively, the rectangular 
chart. The curve corresponding the Cartesian coordinate defined these 
two intersections indicates relative runoff. brief example shows the procedure. 
Let CN, 77; CN, 84; and 3.5 in. (9.89 cm). Using and 3.5 
enter the chart from the top and and enter from the left, relative 
discharge 1.36 interpolated between the curves. The change increases 
discharge 36%. in. (15.24 cm), the relative change due curve number 
would 20%. 

Two limiting lines can identified Fig. The first, vertical line running 
upward from the point coincidence the right, represents relative discharge 
infinity. The second, 45° line running downward and the left from the 
same point coincidence represents relative discharge zero. Points 
and the right this line violate the constraint given with Eq. 


Summary 


The graphical solution for relative discharge provides rapid means for 
estimating runoff changes due watershed modification. The graphical solution 
allows visual appraisal the sensitivity runoff change both modification 
and storm precipitation. Eq. may used for more accurate calculation, but 
should noted that any analysis based curve numbers best 
approximate solution. For problems where actual discharge can measured 
estimated alternative methods, graphical estimation relative change 
provides pragmatic basis for anticipating runoff modification. 
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The following symbols are used this paper: 


SCS runoff curve number; 
storm rainfall, inches (cm); 
direct storm runoff, inches (cm); 
relative runoff; and 
potential maximum retention, inches (cm). 


YOUNG’S INEQUALITY FOR THREE-DIMENSIONAL 
FLows 


INTRODUCTION 


The numerical treatment three-dimensional flows porous media requires 
careful approach regarding the stability and computer memory-time characteris- 
tics the method employed. When extensive information about the flow 
needed, especially for time dependent flows, the computations can quite 
lengthy, and, far the author knows, there are stability criteria general 
validity concerning the numerical procedures used. However, certain cases, 
one may not interested such detailed information, but, some relations 

Center the Public Power Corp., Agion Apostolon St., Athens 811, 
Greece. 

Note.—Discussion open until May 1982. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
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1980. This paper part the Journal the Irrigation and Drainage Division, Proceedings 
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general character concerning the principal flow parameters and consequently 
the use elaborated model, detailed information would not the proper 
choice. Such relation the Young’s inequality which, for steady drainage 
through porous strip contained between two parallel ditches, relates the height 
the water table its center, the precipitation intensity, the ditch spacing, 
and the level permanent water horizon (1). 

This note deals with extension Young’s inequality three-dimensional 
flows. Although the extension formulated with regard transient Darcian 
flow through anisotropic prismatic porous region surrounded peripheral 
ditch (Fig. 1), its main practical significance seems the numerical determina- 
tion bounds the precipitation for steady flows, that the maximum 
height the water table above the impervious boundary will not exceed 


permanent water horizon 


permanent water horizon 
FIG. 1.—Flow Conditions Time, 


given value. After the general form the extended Young’s inequality has 
been presented, specialized steady flows orthotropic regions 
rectangular cross section since for this case one can derive closed type formulas 
for the bounds question. 


The physical process the flow, which stands for the porosity the 
point (x,y,z) the medium while k,, k,, and denote the permeability 
coefficients with respect the directions and respectively, described 
Fig. The other symbols, and appearing Fig. have the meaning 
velocity potential, height water table above the impervious boundary, 


peripheral ditch | 
b x 3 | : 
impervious boundary 
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and precipitation the point (x,y) time, respectively. 

For establishing the extended Young’s inequality the form presented 
this note, additional assumption, regarding k,, and seems necessary 
namely that k,, k,, and are functions and only. Since the form 
this inequality depends the quantity, z)dz, which plays key role 
finding for steady flows bounds associated with maximum height 
the water table, the governing pde z)dz first recorded. 

After slightly modifying the analysis (2) take into account the fact that 
and are not constant and that there also time varying precipitation, 
one can prove that the pde question 


z)dz 
— 


H,, the integral z)dz has the value 


FIG. 2.—Variation versus 


Since the derivation the extended Young’s inequality will based 
geometrical interpretation observation related the physics the 
flow, rough sketch the variation versus while both and are 
kept constant, presented Fig. 

The curve sketched Fig. corresponds the monotonic 
decrease the z-component fluid velocity from h-z condition 
that considered true for the simple drawdown and recharge problems considered 
(2). From Fig. evident that 
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The latest relation first integrated from and then the three 
sides are divided the positive quantity z)dz that finally 


the other hand, since the relation [x,y,h(x,y,t),t] holds, 
follows that 


+ — — 
— 
oy oy z=h oy 0z 


Taking and from Eqs. and respectively, and inserting 
them relation 


— — — = 

which express the conditions satisfied the free surface, not difficult 
prove that 


+ 


The combination Eqs. and leads immediately 


= 
A ah 2 


which constitutes the extended form Young’s inequality three-dimensional 
transient flows through prismatic porous regions surrounded peripheral 
ditch. For steady flows and for the point which corresponds the maximum 
takes the form 
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case one-dimensional steady flow with constant precipitation, through 
isotropic medium contained between two parallel ditches spacing, the 
pde (1) greatly simplified 


+a/2, 


H? 
one transformed 


which constitutes the already known form Young’s inequality (1). More 
generally, for steady flows through media arbitrary cross sections, one can 
solve means the FEM the pde (1) with z)dz the boundary 
and for prescribed functional dependences k,, k,, and and Since 
some cases may not constant, the spatial variation supposed 
have the form f(x,y) with f(x,y) dimensionless function 
intensity. 

After the values z)dz selected points the medium’s cross 
section have been calculated, possible determine approximately the position 
characterized z)dz max which corresponds the maximum height 
the water table well. For f(x,y) and since z)dz 
constant equal the boundary, evident that relation Eq. 
equivalent 


max 


which G(x,y) denotes the solution Eq. which takes the zero value 
the boundary, and the value f(x,y), the point where the maximum 
value G(x,y) has appeared. 

Since the function G(x,y) itself proportional the above relation 
can used establishing bounds the maximum precipitation intensity 
that the height water table will not exceed 

For doubly symmetric problems, i.e., for regions whose geometrical and 
physical properties well precipitation intensity are distributed along two 
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axes symmetry, apparent that appears the center the medium’s 
cross section. For these problems possible apply the FEM the region’s 
quarter only employing simple triangular elements the class (3). 
Example.—Since the main practical use the extended Young’s inequality 
the numerical determination bounds that the height water 
table will not exceed given value, seems worthwhile examine case 
for which closed type formulas which these bounds can found. This example 
corresponds steady drainage through orthotropic region rectangular 
cross section with sides and (see the dotted rectangle Fig. and for 


constant precipitation, i.e., f(x,y) For k,, k,, and constants, 
the pde (1) has the form 


while the solution, G(x,y), satisfying G(x,y) the boundary, 


G(x,y) 


Due the double symmetry the problem, obvious that the maximum 


value G(x,y) appears the center the region 0), that after 
defining 


(1) (2) (3) (4) (6) 
0.100 0.0125 0.400 0.0480 0.700 0.0686 
0.125 0.0156 0.425 0.0504 0.725 0.0695 
0.150 0.0187 0.450 0.0527 0.750 0.0703 
0.175 0.0219 0.475 0.0549 0.775 0.0710 
0.200 0.0250 0.500 0.0569 0.800 0.0716 
0.225 0.0281 0.525 0.0588 0.825 0.0721 
0.250 0.0311 0.550 0.0606 0.850 0.0726 
0.275 0.0341 0.575 0.0623 0.875 0.0729 
0.300 0.0371 0.600 0.0638 0.900 0.0732 
0.325 0.0400 0.625 0.0652 0.925 0.0734 
0.350 0.0427 0.650 0.0664 0.950 0.0736 
0.375 0.0454 0.675 0.0675 0.975 0.0736 
1.000 0.0737 
| 
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(= 


one can set 


are recorded, which were calculated after simple program 
FORTRAN was written. 


Using Eq. for [G(x,y)] 
one takes the form 


max Eq. and noting that the latest 


max 


from which bounds the precipitation intensity can found. 


The main practical significance the extended Young’s inequality seems 
the determination bounds the precipitation intensity for steady flows 
that the maximum height water level will not exceed prescribed value. 
most cases, this made possible after the pde satisfied z)dz 
determined. 

For steady drainage through regions more less cross section 
(e.g. rectangular) and for such form precipitation intensity that formal 


solution the pde (1) possible, closed type formulas for these bounds can 
found. 
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RELATION BETWEEN NUMBER 
AND RUNOFF COEFFICIENT 


INTRODUCTION 


The rational method (1) probably the one most widely used for estimating 
peak discharges for inlet designs small watersheds. The basic input require- 
ments include the drainage area, rainfall intensity for the design return period, 
and runoff coefficient (C). usually applied, the rainfall intensity 
for storm duration equal the time concentration, which dependent 
the land use, the watershed slope, and the hydraulic length the inlet. 
The runoff coefficient most often obtained from table function 
land use, with the table providing range values for each land use. While 
not used widely, tables recommended values exist that allow for variation 
with respect soil type, rainfall intensity, slope, design frequency, and 
surface roughness (3). When range values given for particular land 
use, assumed that the larger values the runoff coefficient should 
used for the larger events (i.e., return periods greater than yr) and large 
slopes; however, the tables often lack sufficient detail exactly what value 
should used any given case. 

The Soil Conservation Service TR-55 graphical method (2,4) another widely 
used technique for estimating peak discharges for small drainage areas. Input 
includes the drainage area, the precipitation depth for the design return 
period, the runoff curve number, and the time concentration. The runoff 
curve number (CN) index runoff potential and function the 
soil type and land use, the soil cover complex ard condition, and the antecedent 
soil moisture. design, the antecedent soil moisture adjustment usually not 
used change the input value the CN. 

general, evident that the two methods, the rational formula and the 
TR-55 graphical method, are based the same input data: drainage area, 
rainfall index, the time concentration, and runoff index. Therefore, may 
interest examine the relationship between the methods. The objective 
this study demonstrate that the two runoff indices (i.e., the curve 
number and the runoff coefficient) are related. 

Dean, Coll. Engrg., Univ. Maryland, College Park, Md. 

Research Asst., Dept. Civ. Engrg., Univ., Maryland, College Park, Md. 
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The two runoff indices, the rational and the SCS CN, are, for numerous 
reasons, important input parameters. First all, and probably foremost, estimated 
peak discharges are very sensitive these indices and any errors estimated 
value will appear error the estimated peak discharge. Secondly, 
comparison between and may reveal differences the rationality 
the two models. Thirdly, both models lack extensive empirical validation. 
The internal validity the rational method has been questioned many. Thus, 
comparison the two methods may identify potential sources bias 
some imprecision the two methods. This would especially useful 
regulatory agency that contemplating change policy that would require 
one method used place the other. Fourth all, this analysis would 
useful when design policies are odds, with the policy one agency 
requiring use one method and another agency requiring designs using the 
other method. The relationship between and then could used adjust 
the (or CN) table that inconsistent design estimates could avoided. 
Fifthly, the relationship would value evaluating policies concerning 
the use different methods for different size watersheds. For example, 
Maryland SCS methods are recommended for most analyses but the rational 
formula permitted for inlet areas less than acre. analysis the relationship 
between and would useful identifying significant discrepancies 
design peak discharges the drainage area point’’ (i.e., the acre 
limit). Sixth all, the two methods use different land use category systems. 
some cases, designer may select one method over another because 
the presentation the land use data; also, the versus relationship may 
helpful refining land use categories eliminate inconsistencies design 
due solely vague land use descriptions. 

The rational formula given 


which the peak discharge, cfs; the dimensionless runoff coefficient; 
the rainfall intensity in./hr for the design frequency and duration 
equal the time concentration; and the drainage area, acre. The 
graphical method given 


which the peak discharge, cfs; the unit peak discharge, 
cfs/sq mile/in. runoff; the drainage area, mile; and the 
runoff volume, inches, computed 


1,000 
which the runoff curve number. The unit peak discharge obtained 
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from graph versus the time concentration, hours (4); this graph 
can represented exponential relationship over small ranges the 
time concentration: 


which and fitting coefficients. Similarly, the rainfall intensity 


- 


100 
CURVE NUMBER (CN) 


FIG. 1.—Relationship between Runoff Coefficient and Curve Number for Baltimore, 
Md. 


the two methods are assumed equally accurate and provide the 
same true value the peak discharge, then the following equality obtained 
from the above equations: 


640 P+0.8S 
Solving for yields 


Since function the CN, then the relation between and can 
obtained for any design frequency, location, and time concentration. Values 
for and are function only the time concentration, while and 
vary with the design frequency, location, and the time concentration. 
Using rainfall intensity-duration-frequency data for Baltimore, Md., and times 


1.0 
0.8 
0.6 
0.4 
4° 
0.2 
0.9 
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concentration two hr, Eq. was solved for return periods (T) 
10, and 100 yr; the resulting values are shown graphically Fig. The 
variation due change time concentration was very small comparison 
with the change return period, and thus the curves Fig. are adequate 
for all times concentration less than hr. The curves Fig. were used 
estimate the runoff coefficient for the range curve numbers across soil 
groups. Table shows the land use categories usually used for the rational 
method and the Soil Conservation Service land use category that most closely 
approximates them. The CNs for soil groups and are given the 
range. These values are used estimate range runoff coefficients from 
Fig. Frequently used values are also shown. For some land uses, the 


TABLE 1.—Comparison Accepted Runoff Coefficients with Values Estimates for 
Curve Number Ranges 


Land Use Values for Range 


Rational Graphical Using Accepted 


(1) (2) (3) (4) (5) (6) (7) 


Business Commercial 
Downtown business 
Neighborhood 

Residential 
Single-family 


0.70-0.95 


Residential 


25% imper- 
vious 
38% imper- 


0.18-0.69 0.32-0.79 


(detached) 
Multi-units 
(attached) 
Residential 
(suburban) 
Industrial 


vious 
65% imper- 
vious 
20% imper- 
vious 
Industrial 


0.43-0.79 


0.52-0.82 


0.53-0.84 


0.14-0.66 


0.61-0.86 


0.66-0.91 


0.28-0.78 


0.73-0.93 


0.60-0.75 


0.25-0.40 


Light 

Heavy 
Parks, ceme- 

teries 


0.50-0.80 
0.60-0.90 
Open space 0.10-0.25 
(good condi- 
tion) 
Open space 
(fair condition) 


0.12-0.71 


Playgrounds 0.25-0.78 0.20-0.35 


two methods agree quite closely runoff potential. For example, the commer- 
cial/business values range from 89-95, which corresponds values 
0.70-0.96 over the three return periods. The accepted values for downtown 
business areas range from 

Similarly, close values result for the industrial land use category. However, 
some land use categories show significant differences. For open space such 
parks, cemeteries, and playgrounds, the range results much wider 
variation values than the recommended values. This may reflect different 
levels impervious area associated with different specific uses within the open 
space category. The residential land use categories also show significant variation. 
This may the result the difficulty associating land use categories. The 


89-95| 0.70-0.87 0.77-0.90 
| | 54-85 

77-92 
51-84 
81-93 

4 = 
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residential land use separation for the rational method too vague that 
one reasonably can expect larger differences runoff potential for single-family 
residential areas than would suggested accepted values (6). Thus, the 
separation would preferable. 


conclusion, the two peak discharge design methods are based similar 
input data and, thus, represent design alternatives; therefore, comparison 
order. The analysis the relationship provides several conclusions: The 
two methods appear give similar peak discharge estimates when one can 
find agreement the land use category separation. For regulatory agencies 
involved specifying the rational formula, appears that consistency design 
estimates could vastly improved better defining the land use categories 
and identifying the other factors affecting variation the runoff coefficient. 

Secondly, peak discharge estimates the rational formula are more sensitive 
values than are estimates the TR-55 graphical method the CN, especially 
for the larger return periods. Increased imprecision that results from this sensitivity 
can controlled better delineation the values within the range 
values the tables. comparison with tables commonly available, the tables 
Rossmiller (3) provide better definition values. The relationship between 
and CN, provided Fig. can used indicate the effects due 
changes soil type, antecedent soil moisture, and land use condition. 

addition the graphical method based the conceptual framework 
the RUNOFF subroutine the SCS TR-20 computer program (5); provides 
reasonable approximation except when the initial abstraction significantly 
greater than about 25% the precipitation. Thus, peak discharges from the 
rational formula also should comparable with those made with the RUNOFF 
subroutine TR-20. 

Fourthly, the values Table estimated from the values, are 
characterized much higher variation than the accepted values. This results 
part because the disparity land use categories. However, many cases, 
the upper value results wide disparity between the estimated and accepted 
values. This partially represents difference the conceptual basis the 
model and would indicate that there would systematic bias peak discharges 
estimated with the two models. When design projects are controlled different 
policies (e.g., watershed size limitations, county state boundaries), this source 
bias should recognized and handled systematic manner. general, 
appears that for residential and open space land use categories the graphical 
method will provide higher peak discharge estimates than the rational method. 
the very least, the relationship can used adjust values for storm 
frequency and soil type; such adjustments may reduce inconsistency estimates. 

also appears that the effect locational differences values could 
assessed evaluating the relationship for locations other than Baltimore; 
however, initial analysis shows that location minor factor, and the results 
for other locations will quite similar. 

And finally, the SCS methods, including the TR-55 graphical method, are 
often criticized because they are based storm distribution. The results 
this analysis indicate that this assumption does not introduce bias into 


4 
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the estimated peak discharge; for those familiar with the way which the 
SCS storm distributions are constructed, the lack bias rational. 
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Soil Conservation Service, Computer Program for Project Formulation: Hydrology, 
Technical Release 20, U.S. Department Agriculture, 1965. 

Whipple, W., Jr., DiLovie, M., and Pytlar, T., Jr., Potential Streams 
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IRRIGATION PRACTICES FOR RECLAIMED 
ALKALI SOILS 


INTRODUCTION 


Alkali soils, about 2,500,000 ha, the Indo-Gangetic alluvial plains are 
characterized impermeable surface soils and hydraulically impeded subsurface 
horizons. This mainly due the high exchangeable sodium (pH 10.4) present 
them. These areas are now being reclaimed farmers large scale 
land levelling and field bunding, applying gypsum (as amendment), and 
cultivating rice and wheat rotation (1,4). 

survey the irrigation practices adopted the farmers the area revealed 
that more than 70% them obtain their irrigation pumping the good quality 
ground water (EC 600 micromhos/cm) and applying average about 
150 and for rice and wheat, respectively, comparison, with 100 and 

Div. Engrg., Central Soil Salinity Research Inst., Karnal-132001, Haryana, 
India. 
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applied the research farm the Central Soil Salinity Research Institute 
(CSSRI), Karnal. Thus, there seems considerable scope for improving 
irrigation efficiencies the farmers’ fields. 


Objectives 
Reclaimed alkali soils require light and more frequent irrigations for wheat 
crop (2). The farmers not provide any flow regulating devices meet this 


FIG. Sloping Borders for Rice Crop 


requirement. This can, therefore, only satisfied minimizing the application 
times through selection proper dimensions (length and width) the border 
strips (Fig. and cutting off the stream when has advanced adequately 
the border. Another possibility adopt sloping border instead 
level one. The sloping border can also adopted rice cultivation 


| 
FIG. 1.—Sloping Borders for Wheat Crop 


402 DECEMBER 1981 


compartmentalization segmentation (Fig. 2). These ideas were evaluated 
series three field experiments CSSRI during 1975-78 with the following 
objectives: (1) Comparison graded and level borders; (2) Determination 
the stream cutoff point; (3) Determination the optimum border dimensions; 
and (4) Adoption sloping borders for rice cultivation. 


Details Experiments 
Experiment No. (1975-76, 1976-77).—Level and graded borders: 


Length, L—90 

Width, 

Slope, S—A level border and border with 0.1% (the prevailing natural slope, 
just smoothed the surface). 

Irrigation cutoff, C,—Four cutoff treatments, viz. when the stream has 
advanced 55%, 70%, 85%, and 100% the total length the border. 

Irrigation depth—The irrigation depth, each irrigation computed 
from the application time and stream discharge. 

Stream discharge—A constant discharge Lps (as obtained normally 
the farmers’ fields) was maintained during each irrigation regulatory 
valve. Parshall flume was installed below the valve for measurement 
discharge. 


Irrigation interval, days.—If were the depth the previous irrigation, 


TABLE Efficiencies and Wheat Crop Yield Graded and Level Borders 
with Different Cut-off 


Graded Border, 
Particulars Level Border Tenths Percent 


Cut-off ratios 
(1) 

Total irrigation water 
applied, centime- 
ters 

Portion the border 
covered each ir- 
rigation, per- 
centage 

Number irrigations 

Application efficiency, 
percentage 

Distribution efficiency, 
percentage 

Crop yield, quintals 
per hectare 

Water use efficiency, 
kilograms per 
hectare per centime- 
ter 


(3) (4) (7) (8) (9) 
values for both years the study. 


Ww 
a 
w 


403 
° 
< 
4 
= 
° 


404 DECEMBER 1981 


the next irrigation will due after interval days, during which the 
cumulative pan evaporation (in cm) equal D/0.7 (8). there any 
rainfall during this period, this amount was added 

Crop—Wheat (S-1553) with recommended fertilizer and agronomic practices 
(1). 

Experiment No. (1976-77, 1977-78).—Border dimensions: 


Irrigation cutoff, C,.—During the first irrigation the stream cutoff when 
has advanced 85% the total length the border. For the subsequent 
irrigations, the corresponding value 70%. 

Other details for this experiment are the same Experiment No. 

Experiment No. (1975, 1976).—Rice sloping borders. The conflicting layout 
requirements rice and wheat rotation are attempted resolved segmenting 
the borders with small and temporary earthen ridges. The segments are intercon- 
nected with small pipes: 


Length, L—93 

Width, W—8 

Slope, 

Average Irrigation D—8 predetermined value from other experiments). 


Experimental Procedures 

The opportunity time different points along the length the border 
determined from the advance and recession curves and the depth cumulative 
infiltration from the infiltration rate curves. The distribution and application 
efficiencies are computed the manner prescribed Israelson and Hansen 


Distribution (%) 


Application efficiency (%) 


Border width (m) 


FIG. 4.—Variation Irrigation Efficiencies and Wheat Crop Yields Different Sized 
Borders 
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(3). Evaluation the effects segmenting the border strips the rice crop 
was done computing the uniformity ponding LS/2nD (6). all 
three experiments, crop yields were recorded study the crop responses. 


Results and Analyses 

Level and Graded Borders (Experiment 1).— Data application and distribution 
efficiencies are presented Table Wheat crop yields and water use efficiencies 
under different cutoff treatments for level and graded borders also are included 
the same table. general, the graded borders are relatively more efficient 
than the level ones. The irrigation data and crop responses also show that 
the graded border with 70% cutoff gives the highest application, distribution, 
and water use efficiencies 84%, 89% and 135.8 kg/ha/cm. The corresponding 
values for the best level border (C, 0.85) the study are 53%, 60% and 
92.1 kg/ha/cm. 

Border Dimensions (Experiment 2).—The values opportunity time selected 
points along the length the border strip, computed from the advance and 


TABLE 2.—Effect Compartmentalizing Graded Borders Rice Yields 


Uniformity Rice Yield, quintals 
ponding per hectare 
Treatment Uc, percent 
(1) (2) 

single compartment. 
border two compart- 

ments. 

compartments. 


recession curves for each border, are shown Fig. The relationships 
are represented nearly horizontal lines, indicating the relative uniformity 
the opportunity time these borders. The variations the mean efficiency 
values and the crop yields (Fig. also confirm this trend. 

Rice Sloping Borders (Experiment 3).—The values the uniformity 
ponding and the rice crop yields are presented Table The results have 
shown that the rice crop yields and values are significantly higher the 
level border strips divided into three compartments (7,) compared 
with 


From these studies possible draw the following conclusions with respect 
the irrigation practices for growing rice and wheat alkali soils: 


layout closed and graded borders can adopted for rice and wheat 
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rotation. The prevailing gentle land slope 0.1-0.2% can immediately adopted 
the border grade. There need attempt for absolute level plots and 
only surface smoothing sufficient. 

Light and frequent irrigations can applied cutting off the irrigation 
when has advanced 70% the border length. 

farmers’ fields, irrigation efficiencies can improved adopting 
graded borders (the prevailing grade 0.1-0.2%) with lengths the range 
50-70 and widths 6-8 These dimensions are particularly suited 
small farms, which are common the region. 

During the rainy season, rice also can grown these graded borders 
dividing them into two three compartments connected asbestos pipes. 
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TECHNICAL NOTES 


DELAYED DRAINAGE STREAM-AQUIFER SYSTEM 
Shlomo 


INTRODUCTION 


recent paper published this journal, Higgins (5) considered the problem 
nonsteady flow from unconfined aquifer into fully penetrating stream 
which the water level fluctuates prescribed fashion. was able 
derive approximate analytical expressions for drawdown any point the 
aquifer assuming that the aquifer material rigid, the water incompressible, 
and the water table moving material boundary. check the validity 
his theoretical expressions, Higgins conducted laboratory experiment with 
the aid sand box. The experimental results were agreement with theory 
for the case where the water level the stream declined linearly with time. 
However, when the water level the stream was lowered instantaneously 
its final position, the experimental drawdown the water table early time 
was higher than that predicted theory; large values time, some 
the experimental water table drawdowns were less than those predicted theory. 
Higgins attributed this lack correspondence between experiment and theory 
failure the part his theory account for unsaturated flow above 
the water table. 

The purpose this note develop alternative theory which does not 
account for unsaturated flow and yet, able explain Higgins’ experimental 
results. The proposed theory similar that Higgins, with two important 
exceptions: (1) The aquifer material assumed elastic and the water 
compressible; and (2) the aquifer taken anisotropic. Whereas Higgins’ 
theory can viewed two-dimensional version the well-flow problem 
reported Boulton (2), the new theory described here derives from the concept 


delayed drainage applied the problem unconfined flow well 
Neuman (6,7). 


Consider unconfined aquifer infinite lateral extent resting imper- 
meable horizontal layer, such that shown schematically Fig. The aquifer 
material uniform but anisotropic with principal permeabilities oriented parallel 
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FIG. 1.—Scheme Stream-Aquifer System 


the coordinates. The aquifer completely penetrated stream infinite 
length. Initially, the system static equilibrium, the water level the 
stream being the same elevation the horizontal water table the aquifer. 
time, the water level the stream suddenly lowered new 
position lying distance, below the original one. result the hydraulic 
gradient generated this drop the water level, water released from storage 
the aquifer due three mechanisms: compaction the aquifer material, 
expansion the water, and gravity drainage the free surface. 


analogy Neuman (6), the problem can expressed mathematically 


a’s 


z 


permeability; specific (elastic) storage; horizontal coordinate; 
vertical coordinate; time, initial saturated thickness; and specific 
yield (drainage porosity). After applying Laplace and infinite sine transforms 
Eqs. 1-6, the resulting one-dimensional problem easily solved. Application 
the inverse sine transform together with Mellin’s inversion formula the 
manner outlined Neuman (6) leads the following solution 


n=0 
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and the terms and are the roots the equations 


. 


(2n 


there gravity drainage (S, the aquifer artesian and 
Note that can also interpreted mean S,b and, thus, 
whereas Eqs. and show that u,(y) for all This 
together with formula from Gradshteyn and Ryzhik (4) and Abramowitz and 
Stegun (1) yields 


which identical with solution developed earlier Stallman (3) for confined 
aquifer. 


For large values one can follow argument identical that Neuman 
(7) show that 


cosh (y) 


which extension Higgins’ (5) equation for the case anisotropic 
aquifer. According Higgins, the solution for isotropic aquifer small 
values time can written (in our notation) 
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form 


Explanation Higgins’ Experimental Results.—It will recalled that small 
values time, the water table drawdowns observed experimentally Higgins 
(5) exceeded those predicted his theory. example, let consider 
his experiment which corresponds 0.33; 0.48 and 0.3. 
From his table find that s,,/c 0.16 and 0.50 
0.030 which s,, drawdown the water table. Eq. due 
Higgins underestimates these values, yielding 0.03 for 
use our own Eq. with 0.2, calculate 0.16 for 
0.017 and 0.30 fort, 0.030, both which are much closer 
the experimental results. The relatively large magnitude required order 
fit Eq. the data may due the presence air bubbles and dissolved 
air the water (7). 

large values time, some the water table drawdowns observed 
Higgins were less than those predicted his theory. This can explained 
simply the fact that and, therefore, Eq. due Higgins will 
generally give smaller values s/c than our Eq. 13. With this, easy 
understand why most Higgins’ data lie above the dashed line, and some 
(including all black dots) lie above the solid line large 
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TECHNICAL NOTES 


SIMPLE EQUATIONS FOR 
EQUIVALENT DEPTH 


INTRODUCTION 


Due its simplicity and relative accuracy compared with other more 
elaborate and hence more complicated methods, the Hooghoudt solution the 
drainage design problem (1) widely used the practicing engineer. addition, 
the notion equivalent depth that Hooghoudt introduced, which makes 
his analysis more important. The equivalent depth has gained considerable 
significance the treatment the transient drainage problem based the 
Dupuit-Forchheimer theory (2,4). 

previous paper the senior writer (3), the solution the so-called 
Hooghoudt equations was generalized the introduction dimensionless 
variables. Application practice for design purposes was facilitated the 
preparation unique nomographic solutions both for the equivalent depth and 
the drain spacing formulas. For calculations computer solving transient 
flow problems, the computation the equivalent depth requires the storage 
the computer memory appreciable size table (maximum 300 values). 

For faster determination the equivalent depth, and avoid storage 
table the computer memory, approximate formulas have been introduced 
(2). Moody’s work was based Hooghoudt’s original work. After the detailed 
reworking the solution Hooghoudt equations (3), was considered 
appropriate proceed with finding simple algebraic equations which approximate 
the equivalent depth more satisfactory way. This the subject the 
present note. 

Previous Relevant the previously mentioned paper (3), the Hoogh- 
oudt equations for computing the equivalent depth and drain spacing were given 
dimensionless form follows: 
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conductivities the pervious soil strata, respectively, above and below the 
plane through the drain axes; maximum height the free-water surface 
above the drain axes level (total loss hydraulic head sustain the ground-water 


flow towards the drains); function the drain radius, given the 
relation 


which the radius the drain; and function the drain spacing 
which given the relations 


2 2 


The value the function (p) can determined easily. The function 
was computed with sufficient accuracy and table has been prepared (3). 
Alternatively, for practical purposes, graph has been constructed (3). The 
table can stored computer memory and thus the determination the 
equivalent depth for solving either Eq. unsteady state drainage problem, 
where the notion equivalent depth utilized, becomes rather easy task. 

Simple Algebraic Equations for Function much simpler approach 
than that previously described for the determination the equivalent depth 
would the use simple algebraic equations approximate accurately the 
function given the complicated relations Eqs. and 3b. Using 
series 100 values for ranging from 0.01-1000, the three curves 


described the following equations were fitted the given data the least 
squares method: For ranging from 0.01-2.00 


For ranging from 2.00-4.60 


4 
"4 


TECHNICAL NOTES 
For greater than 4.60 


Co 


which the prime was used indicate the approximate value 

The error computing the equivalent depth using the approximate values 
the function given Eqs. compared with the values given 

Differentiation Eq. with respect and rearrangement yields the 
relative error due using slightly different values Thus 


The largest relative error occurs for the minimum value which turn 


TABLE 1.—Relative Error Computing Equivalent Depth 


(2) 


occurs for the maximum possible value The latter, designated p,,, 
yields the maximum value for given value Table indicates the 
upper limits the error percentage 

the other hand, estimation the relative error resulting from 
and rearrangement. Thus, one obtains 


confirmed from Tab!e the simple formulas Eqs. approximate 


413 
0.01 —0.08 3.60 
0.10 4.60 0.12 
0.50 4.60 —0.67 
1.00 5.00 —0.40 
1.10 10.00 0.15 
2.00 0.24 50.00 
2.00 100.00 
2.20 0.09 200.00 
2.40 0.20 500.00 0.00 
3.00 0.02 1000.00 0.00 
The largest relative error occurring for 
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the complicated relations Eqs. and closely enough yield error 
the value the equivalent depth less than 1%. 

similar analysis (3) concerning the use the formulas developed 
Moody, the maximum relative error equivalent depth was found 
the order 7%. 


Summary 


Based previous general and dimensionless treatment the Hooghoudt 
equations for computing the equivalent depth and drain spacing, simple algebraic 
relations were developed for easy and sufficiently accurate approximation 
the function Thus, the task computing the equivalent depth for 
solving steady unsteady state drainage problems either desk calculator 
digital computer becomes very simple, while the values computed suffer 
appreciable loss accuracy. 
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The following symbols are used this note (all variables are dimensionless): 


equivalent depth; 

maximum water-table height above drain-axes level; 
hydraulic conductivity the layer above drain-axes level; 
hydraulic conductivity the layer below drain-axes level; 
radius tile drains; 

maximum realistic value for given value 

drain spacing; 

function (Eq. 2); 

function approximating (Eqs. 4). 


K; = 

,(p) = 
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DISCUSSION 


Note.—This paper part the Journal the Irrigation and Drainage Division, 
Proceedings the American Society Civil Engineers, ©ASCE, Vol. 107, No. IR4, 
December, 1981. ISSN 0044-7978 /$01.00. 


415 


16685 IR4 


Discussions 


Discussions may submitted any Proceedings paper technical note published any 
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which have been published ASCE. Discussion paper/technical note open 
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name should indicated below the title (see Discussions herein example) together with 
his ASCE membership grade (if applicable). 

The discusser’s title, company affiliation, and business address should appear the first 
page the manuscript, along with the Proceedings paper number the original 
note, the date and name the Journal which appeared, and the original author’s name. 

Note that the discusser’s identification footnote should follow consecutively from the original 
paper/technical note. the paper/technical note under discussion contained footnote numbers 
and the first Discussion would begin with footnote and subsequent Discussions would 
continue sequence. 

Figures supplied the discusser should designated letters, starting with This also 
applies separately tables and references. referring figure, table, reference that 
appeared the original paper/technical note use the same number used the original. 
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the Publications ASCE for more detailed information preparation and submission 
manuscripts. 


+ 
416 


DISCUSSION 


POLYSTYRENE SPHERES EVAPORATION 


The authors have presented description experiments applying expanded 
polystyrene beads reduce evaporation. When adopting new material and 
methods, indoor and outdoor tests are very useful, not solely prove the 
practical application and its economic advance but for the possibility compari- 
son and for solving theoretical problems, any. this respect, the writer 
would like expand the completion the critical considerations. 

Comparison Evaporation Control Percentages.—According the authors, 
combined cover suppressed evaporation 70.2% and polystyrene 
covers (p. 168). more detailed comparison these figures 
noteworthy, far the contact surface ratios the different cover materials 
are concerned. make more perceptible the evaporation reduction ratios, the 
writer has taken the liberty drawing Fig. which the cross sections 
and the view the used cover material are shown. 

assumed that the contact surface between the polystyrene spheres and 
the water (Fig. can expressed 


This results reduction ratio 42.2%. Thus, the uncovered surface 
0.578 The reduction rate monomolecular cover, according the 
authors’ data, 56.8%. This means, that the case combined cover, 
the evaporation the area covered monomolecular cover 
material will reduced respect the whole area f,, 0.578 
and 56.8 32.7%. Thus, the total reduction, with combined cover, will 
42.2 32.7 74.9%, which approximate the measured value 
reported the authors. The rate somewhat higher than the figure 
70.2%. The deviation most likely arises due distance between spheres and 
varied sphere diameters. 

Influence Dimensions and Specific Gravity Sphere.—It was reported that 
beads varied mm-8 diameter. the case polystyrene spheres, 
the rate evaporation influenced some extent the different diameters 


and shell thicknesses the spheres and the value specific gravity 
that material. Namely: 


the diameter, the sphere larger, the water displacement due 
net weight should have relatively smaller section diameter which will 
decrease the suppression ratio. 


1980, Tariq Mahmoud and Natiq Bashi (Proc. Paper 15689). 
“Tech. Advisor VITUKI, Research Centre for Water Resources, Budapest, Hungary. 
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the specific gravity higher the sphere will submerge deeper the 
water, consequently, will larger and the reduction ratio will increase. 


Different Temperature Elevations Water.—The tested polystyrene spheres 
diameter byproduct whose material has considerable insulation 
features. This one the reasons for the temperature increase. the other 
hand, one should take into consideration the fact that, during evaporation, the 


Polystyrene spheres 


Contact surface 


Section a—a 


Circle section f of 
contact surface 


Top view 


FIG. 8.—Cross Sections and View Used Cover Material 


vapor molecules leaving the water surface transfer certain heat, the 
called i.e. 


which the water temperature, degrees Celsius. The characteristic 
line water temperature running down the water surface has been studied 
detail Kristina Katsaros, al. (13). The writer has adopted these 
experimental characteristic diagrams the case water covered polystyrene 
spheres. The full line, marked Fig. shows the water temperature diagram, 


aa \ / \ 
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plotted against the depth uncovered water The dotted line, marked 
Fig. shows the same diagram for water covered with polystyrene spheres. 

Between the two lines are the kind differences. one respect the 
evaporation rate reduced the cover material, the heat loss water 
lower, consequently the surface temperature will higher 


the other hand, the insulation features the polystyrene spheres 
used, the heat convection reduced which has been resulted higher (1.6° 


Polystyrene spheres 


FIG. 9.—Water Temperature versus Depth Uncovered Water 


6.5° 4.2° water temperature. This water temperature elevation marked 
Fig. generally with 


15. Katsaros, B., Liu, T., Businger, A., and Tillman, E., Transport 
and Thermal Structure the Interfacial Boundary Layer Measured Open Tank 
Water Turbulent Free Journal Fluid Mechanics, Vol. 83, Part 
London, England, Nov., 1977, pp. 311-335. 
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